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ABSTRACT 

This  thesis  reports  on  a  simulation  study  of  parametric  and  nonparametric 
procedures  for  obtaining  confidence  intervals  for  the  logarithm  of  the  probability  a 
serru-markov  process  enters  a  particular  state  before  a  fixed  time  t.  Three  estimators 
and  confidence  interval  procedures  are  proposed  and  compared.  The  different 
estimators  use  different  amounts  of  information  about  the  process.  The  maximum 
likelihood  estimator  and  its  normal  confidence  interval  procedure  uses  the  most;  the 
estimator  based  on  the  empirical  distribution  function  of  the  observed  first  passage 
times  uses  the  least.  An  estimator  based  on  an  exponential  approximation  to  the 
survivor  function  of  the  first  passage  time  uses  an  intermediate  amount  of  information; 
confidence  intervals  for  the  last  estimator  are  obtained  using  jackknife  and  bootstrap 
procedures.  The  maximum  likelihood  procedure  is  the  most  efficient  if  the  underlying 
model  is  correct.  If  the  model  is  not  correct  the  empirical  survivor  function  estimator 
appears  to  be  best  for  small  times  and  the  estimator  based  on  the  exponential 
approximation  best  for  large  times. 
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I.  INTRODUCTION 

A.  OBJECTIVES 

Finite  state  space  semi-Markov  models  find  application  in  a  variety  of  areas  such 
as  queueing  theory,  reliability,  and  clinical  trials.  The  application  of  these  models  often 
centers  on  the  distribution  of  a  first-passage  time  to  a  state  or  a  set  of  states 
representing  for  example  the  lifetime  of  a  system  or  marking  the  end  of  a  busy  period 
for  a  server.  Suppose  that  the  observations  of  the  path  of  the  semi-Markov  process  are 
all  that  is  known  about  the  process. 

This  thesis  reports  the  results  of  a  simulation  experiment  to  compare  various 
parametric  and  nonparametric  estimators  of  the  natural  logarithm  of  the  probability  a 
semi-Markov  process  does  not  enter  a  particular  state  before  a  given  fixed  time  t.  In 
what  follows  we  will  use  "In  probability"  for  "  natural  logarithm  of  the  probability  ". 
The  specific  semi-M-arkov  model  and  estimators  considered  are  given  in  Chapter  II. 
Chapter  III  contains  the  details  of  the  simulation  experiment  and  results.  Conclusions 
from  the  study  are  given  in  Chapter  IV. 

B.  SCOPE  OF  THE  THESIS 

The  purpose  of  this  thesis  is  to  use  simulation  to  compare  estimators  and  their 
confidence  intervals  for  the  In  probability  a  particular  semi-Markov  process  does  not 
enter  a  particular  state  before  a  given  fixed  time  t.  The  particular  semi- Markov  model 
and  estimators  considered  are  given  in  Chapter  II.  A  simulation  study  comparing  bias 
and  standard  errors  for  these  estimators  was  reported  in  Gallagher  (1986)  [Ref.  1].  In 
this  thesis,  we  are  primarily  interested  in  comparing  confidence  interval  procedures. 

An  estimation  procedure  which  uses  the  least  information  about  the  semi- 
Markov  process  uses  the  empirical  distribution  of  the  observed  first  passage  times.  An 
estimation  procedure  which  uses  the  most  information  is  to  assume  a  parametric  form 
for  the  sojourn  time  distribution  in  each  state  and  a  transition  matrix  to  describe  the 
transition  between  states;  the  parameters  of  the  sojourn  time  distribution  and 
transition  probabilities  can  be  estimated  using  maximum  likelihood.  An  estimation 
procedure  requiring  less  information  uses  nonparametric  estimators  of  the  sojourn  time 
distributions  and  the  maximum  likelihood  estimators  for  the  transition  probabilities. 
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In  many  cases,  parametric  assumptions  concerning  the  sojourn  time  distribution 
are  difficult  to  justify.  It  was  demonstrated  in  Park  (1986)  [Ref.  2]  that  incorrect 
parametric  assumptions  may  lead  to  very  biased  estimators.  Hence,  a  nonparametric 
estimation  procedure  may  be  preferred  to  a  parametric  one  when  actual  data  is  used. 
However,  the  nonparametric  procedure  can  be  expected  to  be  less  efficient  than  a 
parametric  one  provided  the  parametric  model  is  correct. 

The  thesis  is  organized  as  follows.  In  Chapter  II,  the  nature  of  the  problem  is 
described  and  several  parametric  and  nonparametric  confidence  interval  estimation 
procedure  are  introduced.  In  Chapter  III,  the  simulation  experiment  is  described  and 
results  are  given.  Conclusions  drawn  from  the  simulation  study  are  given  in  Chapter 
IV. 
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II.  NATURE  OF  THE  PROBLEM 

A.       DESCRIPTION  OF  PROBLEM 

The  semi-Markov  process  model  used  in  the  simulations  to  compare  the 
estimators  is  as  follows. 

Suppose  we  observe  N  individuals.  Let  Xt  (i)  be  the  state  of  the  i-th  individual 
at  time  t.  We  will  assume  [Xt  (i),  t^  0}  1=1,  2,....,  N,  are  independent  semi- 
Markov  processes  with  three  states  (  0,  1,  2  }  having  the  same  distribution  as  [X*  .  t2: 
0}.  All  individuals  start  at  t  =  0  in  state  1,  (XQ  =  1).  An  individual  stays  in  state  1  for 
a  random  length  of  time  having  distribution  function  Fj.  Upon  leaving  state  1.  the 
process  transitions  to  state  0  with  probability  9  and  to  state  2  with  probability  1-9  . 
If  the  process  transitions  to  state  2.  it  spends  a  random  length  of  time  there  having 
distribution  function  F0.  From  state  2,  the  process  transitions  to  state  I  with 
probability  I.   State  0  is  an  absorbing  state.   Once  in  it,  the  individual  never  leaves. 

For  all  individuals,  the  entire  path  of  transitions  and  sojourn  times  are  observed 
until  the  time  of  absorption  in  state  0  (Fig.  2.1). 

Let 

D  =  inf{t>  0,  Xt  =0}, 

be  the  entrance  time  to  state  0  (or  the  time  of  death). 

The  problem  is  to  estimate  the  logarithm  of  the  first  passage  time  survivor 
function  P  (D>t}  for  fixed  time  t  from  the  data  obtained  by  observing  the  N 
individuals. 


B.       THEORETICAL  DEVELOPMENTS 

Two  estimators  for  P^D>t]  will  be  described  in  this  section.  They  are  the 
maximum  likelihood  estimator  and  the  asymptotic  renewal  estimator  from  a  paper  by 
P.A.  Jacobs.    [Ref.  1:  p.2-6] 
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Figure  2.1     Transition  State  Step  (Sojourn  Time). 

1.  Maximum  Likelihood  Estimate  for  Continuous  Time  Markov  Chain 

In  this  subsection,  the  maximum  likelihood  estimator  for  P(D>  t)  will  be  given 
for  the  special  case  in  which  the  sojourn  time  in  state  i  is  exponentially  distributed 
with  mean  l,p-  (i=  1,  2). 

Let  R  be  the  number  of  transitions  from  1  to  2  for  one  individual.  The  log 
likelihood  function  for  the  individual  is 


L  =  R  In  (1-0)  +  In  8  +  R  In  p- 


(eqn  2.1) 


+  (  1  +  R )  In  p1  -  p]Tl  -  p2T2 


where  T-   (i=  1,2)   is  the  total  time  spent  in  state  i  before  death 
The  maximum  likelihood  estimators  are 


A 

e 


A 
Pi 


1 


1  +  R 

1  +  R 


(eqn  2.2) 
(eqn  2.3) 
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A  R 

p2  =  —  (eqn  2.4) 

[2 

Further  R  is  geometric   {0}     so   E{R}    =    (1   -  0)/6.     Since  the   N  individuals  are 
independent,  the  maximum  likelihood  estimators  using  the  data  for  all  N  individuals 
are 


a  N 

G  =  —  (eqn  2.5) 

N  +  R 


a        N  +  R 

?l  =  — ; (eqn  2.6) 


a  R 

P2  =  t-  (eqn  2.7) 


where  R  is  the  total  number  of  transitions  from  1  to  2  for  all  N  individuals  and  T.  is 

i 

the  total  time  spent  in  state  i  for  the  N  individuals.    For  the  estimators  based  on  data 
for  all  N  individuals. 


Var  {$}  =    I  (0)     [  (eqn  2.8) 


Var  {p*}  =  I  (Pl)     l  (eqn  2.9) 


var  (p2}  =  I(p2)     l  (eqn  2.10) 

where 


1(0)  =  N  E  {-1-  +  -^j-}.  (eqn  2.11) 

«  (i-«r 
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I(Pl)  -  N  E  {  \R    }  (eqn  2.12) 

Pi 


I(p2)  =  NE{— j-J  (eqn  2.13) 


with  N  being  the  number  of  individuals. 

Let  D  be  the  time  of  entrance  to  state  0  for  an  individual.    Fix  t  and  put 

S  =  S(t)  =  P  {D>t}; 

then 

^+  P->  .  «        ^i +  P^        Op, 

S(t)  =  {      -,      2  EXP  (k2  i)  -  EX?(X{  t)     l,      -      )  P-  (eqn  2.14) 

.  .  K-y  A.  A .  ■  A.» 

where  Ap  A-,  are  the  roots  of  the  equation 

0P1P2  +  y  (PL  +  P2)  +  y2  =  0  .  (eqn  2.15) 

A  parametric  estimator  for  S(t),  the  survival  function,  is 


A  A  A  AAAAA 

A  k^  +    p.,  K.  +    p^      0p, 

PM  (D>t}=  {    -   a    -    EXP(M  -  EXP^t)    \  Vl  }^  (eqn  2.16) 

A.^  A,  A ,  "A.» 

where  A,  and  A-,  are  roots  of  the  equation 


AAA  A  A  , 

0P,P->  +  y  (P,  +  P-0  +  y~  =  0  .  (eqn  2.17) 


Since  the  maximum  likelihood  estimators  are  uncorrected,  the  asvmptotic 

A  A 

variance  of  P^  [D>  t) ,  which  will  be  denoted  by  Var^j,   is  approximately 

VarM  !PM  !D>t}|0.p1.p2]  (eqn  2.18) 
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=  Var  {6}  (BS:dQ).2  +  Var  {p{}  KdSidp^2  +  Var  ip2V   (<3S  cp,)2 


where 


3S   3  ,  gS   5A1  +  35   3A2 
30   0   3a   30    oXp  39 


(eqn  2.19) 


3S    S    33 


3X-j 


(eqn  2.20) 


jjp1    p.j    oA1  ap 


OA,   OP. 


ep 
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The  asymptotic  variance  of  In  P^  {D>t}  =  Var^  /  (Py  (D>t})ib. 
2.  Asymptotic  Renewal  Estimator 

This  subsection  describes  an  estimation  procedure  for  P  (D>  t),  which  is  based 
on  an  exponential  approximation  to  the  survivor  function  P{D>t}.  Details  of  the 
approximation  which  is  obtained  by  asymptotic  renewal  theoretic  results  is  given  in 
Jacobs  [Ref.  3].   The  approximation  improves  as  t  becomes  larger. 

A 

Let   9  =  N  ,  (N  +  R  )    be  the  maximum  likelihood  estimator  of  0  (equation 

A 

2.5).    Let    (Pj(a)  be  the  empirical  transform  of  the  sojourn  time  in  state  i  ;  that  is,  if 
SA(i),  S2(i), S>yf  (i)  denote  all  the  sojourn  times  in  state  i  for  the  N  individuals,  then 


M. 

A 

q>i  (a)  =  (1/Mj)  I  EXP  (a  x  Sk  (i) ).  (eqn  2.27) 

k=l 
The  asymptotic  renewal  estimator  of  the  survival  distribution  [Ref.  3]  is 


PA{D>t}  =  (b/|l)  EXP(-K  x  t).  (eqn  2.28) 

A 

where   K.  is  the  solution  to  the  equation 


f(a)  =  (1  -$)<p  (a)<p,  (a)-  1  =  0.  (eqn  2.29) 


Mi  M^ 
a      (1-9)      -1-2  a 

H  =  v   v  [Sk(l)  +  S:(2)}  EXP(K(Sv(l)  +  S-(2))}.  (eqn  2.30) 

1         Tt=li  =  l 


q)1(K)=  (1/Mj)  VEXP(KX  Sk(l)).  (eqn  2.31) 

k-1 


and 


A  A    A 

b  =  (0; 


(8/K)  <pj(K).  (eqn  2.32) 
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Note  that  when  a  =  0,  the  left  hand  side  of  the  equation  (2.29)  is  -  0<  0.   As 
a  increases,  the  left  hand  side  of  the  equation  increases  to  oo,  thus  there  is  a  unique 
solution  k.   The  solution  may  be  found  by  numerical  search.   One  possible  numerical 
search  procedure  is  the  golden  section  search. method. 
3.  Empirical  Distribution  estimate 

Let  dp  d-,, ,  d^-  denote  the  observed  times  of  absorption  in  state  0  for  the  N 

individuals.   A  binomial  estimator  for  the  survivor  function  is 


N 
PB{D>t}  =  (1/N)  J  l(t>00)  (d.)  (eqn2.33) 

k=l 


where      l(t?oo)  (dj)  =  *   ifd.>t, 

.    .  0   otherwise. 

C.       CONFIDENCE  INTERVAL  PROCEDURE 

Suppose  X  is  a  random  variable  whose  probability  law  depends  on  an  unknown 
parameter  0.  Given  a  random  sample  of  X:  Xj,  x->,.---,  x  the  two  statistics  lower  (L) 
and  upper  (U)  form  a  100(1  -  a)%  confidence  interval  for  0  if 

P{L  <  0  <  U}  >  1  -  a. 

Procedures  to  obtain  confidence  intervals  for  the  point  estimates  for  In  P{D>t} 
are  given  below. 

1.  Maximum  Likelihood  Estimator 

The  maximum  likelihood  estimator  is  asymptotically  normal  as  the  number  of 
individuals  N  becomes  large,  [Ref.  4].  For  a  fixed  t,  the  maximum  likelihood 
confidence  interval  for  In  P{D>  t}  is 


{  L  ,  U  }  =  In  PM{D>t)  ±zi.aj2  (VVarM)  /  (PM  {D<t}).  (eqn  2.34) 

A 

where  Var^  is  obtained  by  using  the  maximum  likelihood  estimator  of  0,  pv  and  p2  in 
equation  (2.11)  -  (2.26).   Confidence  limits  that  are  larger  than  0  are  set  equal  to  0. 
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2.  Empirical  Distribution  of  First  Passage  Time  (  Binomial  and  Normal  Approx. 
Confidence  Interval   Methods  ) 

a.  Binomial  Confidence  Interval 

Since  the  N  individuals  are  assumed  independent,  the  estimator  Pg{D>t) 
of  equation  (2.33)  has  a  binomial  distribution  with  N  trials  and  probability  of  success 
P(D>  t}.  The  IMSL  routine  BELBIN  was  used  to  obtain  binomial  confidence  intervals 
for  P{D>t}.  For  a  description  of  the  procedure  see  [Ref.  5:  p. 390  -  391].  The  binomial 
confidence  interval  for  In  P(D>t]  is  obtained  bv  taking  loaanthms  of  the  upper  and 
lower  confidence  limits.  If  the  lower  confidence  limits  for  Pg  (D>  t}  is  equal  to  0.  it  is 
set  equal  to  0.0001. 

b.  Approximate  Normal  Confidence  Interval. 

If  the  number  of  individuals  N  is  large,  the  binomial  confidence  interval 
procedure  can  be  approximated  by  a  normal  confidence  interval  procedure  as  follows, 
[Ref.  6:  p.99  -  100  .  p.954  -  955].   The  interval  is 


(  L  ,  U  }  =  PB(D>t}  ^Va^VVarg  (eqn  2.35) 

where 


VarB  =  (  l/N){PB(D>t}x(  [.   PB 


This  interval  {  L  ,  U  }  is   an  approximate  100(1  -  a)%  confidence  interval  for  P[D>t}. 
The  interval  for  In  P(D>  t}  is  {  In  L  ,  In  U  }.  If  either  L  or  U  is  less  than  0,  it  is  set 
equal  to  0.0001.  If  either  L  or  U  is  greater  than  1,  it  is  set  equal  to  1. 
3.  Asymptotic  Renewal  Procedure 

In  this  subsection,  two  nonparametric  methods  are  described  for  obtaining 

A 

confidence  interval  using  the  estimator  P^  ;D>t}.    One  is  the  jackknife.  The  other  is 
the  bootstrap.  They  are  described  below. 
a.   Jackknife  Estimation  Method 

The  jackknife  was  first  introduced  by  Quenouille  (1956)  [Ref.  7],  for  the 
purpose  of  reducing  the  estimate  bias,  and  the  procedure  was  later  utilized  by  Tukey 
(1958)  [Ref.  8],  to  develop  a  general  method  for  obtaining  approximace  confidence 
intervals. 


The  basic  idea  of  the  jackknife  estimation  method  is  to  assess  the  effect  of  each 
of  the  groups  into  which  the  data  have  been  divided,  not  by  the  results  for  that 
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jroup  alone,  but  rather  through  the  effect  upon  the  bodv  of  data  that  results 
rom  omitting  that  eroup.    The  two  bases  of  the 
lesired  captation  for  all  the  data,  and  then,  aft 
,ve  make  the  calculations  for  each  of  the  slightly 
by  leaving  out  just  one  of  the  groups.   [Ref.l?] 

The  jackknife  procedure  is  as  follows: 

A 

(1)  Let  Yan  (t)   be  the  estimate    P^(D>  t}  computed  using  all  the  data. 

(2)  Let  Y:  (t)  be  the  computed  statistic  using  that  data  which  omits  the  j 

subgroup  where  j   =    1,  2 ,  k.    In  this  thesis,  for  the  case  in  which  the 

number  of  individuals  N  equals  20    the  j       subgroup  consists  of  all  data 

f  h 

corresponding  to  the  j  individual.  For  the  case  in  which  the  number  of 
individuals  N  equals  50  the  first  subgroup  consists  of  all  data  corresponding  to 
the  first  5  individuals,  the  second  subgroup  of  the  second  5  individuals  etc. 
Some  cases  for  N  =  50  were  run  with  the  j  subgroup  consisting  of  the  j 
individual.  The  resulting  confidence  intervals  differed  little  from  .those  obtained 
by  leaving  out  5  individuals  at  time.  Computational  considerations  lead  us  to 
use  the  fewer  number  of  subgroups  in  this  case. 

(3)  Define  the  j1"-  pseudo-value  by 

Y*j(t)  =  k  In  Yall(t)  -  (k-1)  In  Y:(t).  (eqn  2.36) 

(4)  The  jackknife  estimator  Y*(t)  is 


Y*(t)  =  (1/k)  {Y*    (t)  +Y*    (t)+ 4-Y*    (t)}.  (eqn  2.37) 

1  2  k 

(5)      The  jackknife  estimator  of  the  variance  of  Y*  (t)  is 


S*2(t)  =   [l/(k  (k-1))}  £  (Y*j(t)  -  Y*(t))2.  (eqn  2.38) 

j=l 


Tukey(1958)  proposes  that  in  a  wide  variety  of  problems  the  k  estimated 
pseudo-values  can  be  treated  as  approximately  independent  and  identically 
distributed   random  variables   [Ref.  8],    to    obtain   the    following   confidence 
interval  procedure. 
(6)      The  jackknife  confidence  interval  is  computed  as  follows. 
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{     L  ,  U  }  -  Y*(t)  ±  tj.a/2  VS*2(t)  (eqn  2.39) 

where  t,  «/,  is  the  upper  1-  a/2  critical  point  of  the  t-distribution  with  k-1 
'    degrees  of  freedom.     The  confidence  interval  given  by  equation  (2.39)  is  a 

function  of  the  estimated  variance.  If  either  confidence  limit  is  greater  than  0, 

it  is  set  equal  to  0. 

b.   Bootstrap  Estimation  Method 

Efron(1979)  introduced  the  bootstrap  method  for  estimating  the 
distribution  of  a  statistic  computed  from  observations,  [Ref.  10].  In  this  thesis,  the 
bootstrap  was  implemented  as  follows. 

Suppose  data  are  gathered  for  N  individuals.  Let  Rp(n)  be  the  number  of 
transitions  from  state  1  to  state  2  for  individual  n.  Let  [S^(i)}  be  the  collection  of  all 
sojourn  times  in  state  i  for  all  individuals.  A  bootstrap  replication  is  generated  as 
follows:  To  generate  data  for  one  individual,  one  observation  is 'drawn  at  random  with 
replacement  from  (Rp(n)};  call  the  observation  r,,;  r,,"*"  1  observations  are  drawn  at 
random  with  replacement  from  the  collection  of  state  1  sojourn  times  (Sj.(l)};  rp 
observations  are  drawn  at  random  with  replacement  from  the  collection  of  sojourn 
times  in  state  2,  (Si,(2)}.  This  procedure  is  replicated  N  times  to  generate  bootstrap 
data  for  N  individuals.  The  estimator  In  P^  [D>t}  is  computed  using  the  generated 
data.    This  completes  one  bootstrap  replication;  B  bootstrap  replications  are  done. 

A 

The  B  estimates  of  In  PA   {D>t}are  ordered.    A  100x(l-a)%  confidence  interval  is 

n  A 

constructed  using  the  a/2  and  1  -  a/2  quantiles  of  the  bootstrap  estimates  of  In  P^ 
{D>t}.  If  either  confidence  limit  is  larger  than  0,  it  is  set  equal  to  0.  In  the 
simulations  the  number  of  bootstrap  replications  is  B  =  100. 
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III.  ANALYSIS  OF  THE  CONFIDENCE  INTERVAL  PROCEDURES 

A.       DESCRIPTION  OF  SIMULATION 

A  Fortran  program  is  written  to  simulate  the  semi-Markov  process  and  compute 
the  confidence  interval.  All  simulations  are  carried  out  on  an  IBM  3033AP  Computer 
at  the  Naval  Postgraduate  School  using  the  LLRANDOM  II  Random  number 
generating  package  [Ref.  13]. 

The  data  for  the  simulated  experiments  are  generated  as  follows: 

(1)  An  individual  starts  in  state  1  at  time  0. 

(2)  A  random  number  with  distribution  F1  is  generated  for  the  sojourn  time  in 
state  1. 

(3)  A  uniform  random  number  is  generated. 

If  it  less  than  9.  the  process  transitions  to  state  0  and  the  data  for  one  individual  is 
complete.  If  it  is  greater  than  9.  the  process  transitions  to  state  2.  A  random  number 
having  distribution  function  F^  is  generated.  The  procedure  then  returns  to  step  2. 
Data  are  generated  for  N  individuals  and  are  collected  as  follows: 

(1)  The  passage  times  to  state  0  for  each  individual; 

(2)  The  sojourn  times  in  states  1  and  2  for  each  individual: 

(3)  The  number  of  transitions  from  state  1  to  state  2  for  each  individual. 

For  each  replication  of  the  simulation  the  estimates  and  confidence  interval  in 
Chapter  II  are  computed  for  In  P{D>t}.  Each  simulation  is  replicated  300  times.  The 
true  In  P(D>t}  is  computed.  For  each  procedure,  the  number  of  confidence  intervals 
covering  the  true  value  is  recorded. 

The  number  of  individuals  that  are  too  low  (  true  In  P(D>t}>LT;  the  upper 
endpoint  of  the  interval  )  and  too  high  (  true  In  P(D>  t}  <  L;  the  lower  endpoint  of  the 
interval  )  are  also  recorded.  The  average  length  of  the  confidence  interval  is  computed 
as  well  as  the  standard  deviation  of  the  lengths.  The  standard  deviation  is  computed  by 
subtracting  the  mean  length  from  each  length,  squaring  the  results,  summing  over  the 
300  lengths  and  finally  divided  by   299. 
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B.       SIMULATION  RESULTS  FOR  THE  EXPONENTIAL  MODEL. 

In  this  section  results  will  be  reported  for  a  simulation  experiment  in  which  the 
sojourn  time  distributions  in  both  states  are  exponential;  Fj  is  exponential  with  mean 
1/pj  and  F2  is  exponential  with  mean  1/  p2.  Some  true  values  for  P(D>  t}  for  this  case 
can  be  found  in  Appendix  A. 

For  each  replication  nominal  80%,  and  90  %  confidence  intervals  for  In  P{D>  t} 
for  various  values  for  t  are  computed  using  each  procedure  of  Chapter  II.  The  times 
considered  are  t  =0.5,  1.0.  1.5,  2.0,  3.0,  4.0. 

The  confidence  interval  procedures  for  In  P[D>t}  are  the  binomial  confidence 
interval  (=  BIN  ),  and  it's  approximating  normal  confidence  interval  (=  NOR  ),  for 
the  In  fraction  of  individuals  who  have  not  entered  state  0  by  time  t:  the  maximum 
likelihood  confidence  interval  (=  MLE  ),  jackknife  confidence  interval  for  the 
asymptotic  renewal  estimator  (=  J.K  ),  and  bootstrap  confidence  interval  for  the 
asymptotic  renewal  estimator  (  =  •  B.T  ).■  For  each  procedure  the  number  of  intervals 
covering  the  true  value  of  In  P  [D>t}  is  recorded  as  well  as  the  number  of  intervals 
that  are  too  high  or  too  low.  These  results  are  reported  in  Tables  1,  3,  5  and  7.  Next  to 
each  coverage  count  are  given  the  corresponding  coverage  proportion  in  parenthesis. 
Below  the  coverage  numbers  is  given  a  95%  confidence  intervals  for  the  coverage  rates. 
This  interval  is  computed  as  follows.  If  P  is  the  proportion  of  (l-a)%  intervals  that 
cover  the  true  value  then  a  95%  confidence  interval  for  the  coverage  rate  is 


{  CL  ,  CU  }  =  P  ±  1.96V  P  (1-  P  )/300. 

If  a  (l-a)%  confidence  interval  procedure  is  performing  well,  then  this  interval  should 
cover  about  (l-a)%  of  the  time.  In  addition,  if  a  80%  (respectively  90%)  confidence 
interval  procedure  is  working  well,  then  out  of  300  replications  between  226 
(respectively  260)  and  254  (respectively  280)  confidence  intervals  should  cover  the  true 
value.  Simulations  for  which  the  numbers  of  intervals  that  cover  are  outside  of  these 
bounds  are  given  in  bold  type.  In  the  first  column  of  Tables  1,  3,  5  and  7  the  true 
?{D>  t}  is  given  in  parentheses. 

In  Tables  2,  4,  6  and  8  are  recorded  the  average  length  of  the  confidence  intervals 
for  In  P{D>t}.  The  estimated  standard  deviation  of  the  length  is  below  the  average 
length  in  parenthesis.  If  a  procedure  is  performing  well,  it  should  not  only  have  the 
correct  coverage  rate  but  also  a  small  average  length.  The  simulation  results  recorded 
in  Tables  1  -  4  are  for  a  simulation  whose  parameter  values  are  p.  =  1.0,  p0=  10.0  and 
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8  =  0.5.    The  number  of  individuals  N  is  set  at  20  and  50,  representing  a  low  and 
moderate  number  of  individuals. 

Coverage  results  for  N  =  50  individuals  are  presented  in  Table  1.  The  binomial 
confidence  interval  procedure  tends  to  overcover.  The  maximum  likelihood  procedure 
has  the  right  coverage  for  all  cases.  The  two  confidence  interval  procedures  using  the 
asymptotic  renewal  estimators  undercover  for  t<  1.0  and  have  the  right  coverage  for 
t>  1.5.  When  a  jackknife  confidence  interval  does  not  cover  it  is  often  because  it  is  too 
high  (true  value  <  L  ). 

Sample  means  and  variance  of  the  confidence  interval  lengths,  for  N  =  50 
individuals,  are  reported  in  Table  2.  The  average  lengths  for  the  maximum  likelihood 
and  bootstrap  are  very  close  for  t  larger  than  1.0.  The  bootstrap  and  jackknife 
procedures  have  larger  average  confidence  interval  lengths  than  those  for  the  maximum 
likelihood  estimator  but  smaller  than  for  the  binomial  procedure. 

Results  for  a  simulation  with  N  =  20  individuals  are  given  in  Tables  3  and  4. 
They  are  similiar  to  those  in  Tables  1  and  2.  However,  the  average  length  of  the 
intervals  is  larger  in  Table  4  than  in  Table  2  retlecting  the  smaller  number  of 
individuals. 

Results  of  a  simulation  with  the  same  parameters  as  the  first  but  with  the  9  = 
2/3  are  given  in  Tables  5  -  3.  Results  for  N  =  50  individuals  are  given  in  Tables  5  and 
6.  The  coverage  results  are  in  Table  5.  The  binomial  procedure  has  better  coverage 
than  in  Table  1.  The  procedure  based  on  the  asymptotic  renewal  estimators  have  the 
correct  coverage  only  for  t  >  2.0.  The  average  length  of  the  intervals  in  Table  6  are  in 
general  longer  than  the  lengths  in  Table  2. 

Results  for  a  simulation  with  9=  2/3  and  N  =  20  individuals  are  given  in  Table 
7  and  8.  The  binomial  confidence  interval  tends  to  overcover  when  the  true  probability 
is  greater  than  0.5241  or  less  than  0.2753.  Once  again  the  confidence  intervals  based 
on  the  asymptotic  renewal  procedure  undercover  for  t  <  2.0.  The  maximum  likelihood 
confidence  intervals  give  the  correct  coverage. 
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TABLE  1 

COVERAGE  RATIO  (EXP  MODEL) 
N  =  50,  8  =  0.5,  Pj  =  1.0,  p2=  10.0 


time  C.I  (%) 


coverage 


BIN 


NOR 


MLE 


J.K 


B.T 


0.5      80  % 


(.7866) 


too  hieh 
cover 

too  low 
int.val 

too  hieh 
cover 

too  low 
int.val 


5(.02) 
288(.96) 

7(.02) 
(.94.-98) 


38(.13) 
239(.SO) 

23(.08) 
(.75,.84) 

17(.06. 
276(.92 

7(.02V 
(.89,.95) 


245 


(.77,.86) 


14(.05^ 
274  .91) 

12  .04) 
(.88  .95) 


119(.40) 
167(.56) 

14(.05 
(.50,.61 

83(.28) 

209(.70) 

8(.03y 

(.64.-75) 


15H.50) 
149(.50) 

0(.0O) 

(.44,.53) 

110(.37) 
190(.63) 
Ojf.OO) 
(.5S..69) 


1.0      so  % 


(.6203) 


90  % 


too  high 
cover 

too  low 
int.val 

too  hieh 
cover 

too  low 
int.val 


28(.09) 

245C.S2 

27(.09 

(,77;.S6) 

8(.03) 
282(.94) 
10(.03) 

(.91..97) 


28(.09) 

245  .82) 

27  .09) 

(.77..S6) 

12(.04) 

275(.92) 
■  13(.04) 
(.S9..95) 


27(.09) 
245(.S2 

28  .09) 
(.77;.86) 

14(.05) 
274  .91 

12  .04 
(.88.951 


70(.23) 

216(.72) 

14(.05) 

(.67;.77) 

45(.15) 
247(.82) 

8(.03J 

(.78,.87) 


80(.27) 
213(.71) 

(.66,.76) 


45(.15) 
250(.83) 

"  02) 

88) 


M 


1.5      80  % 


(.4891) 


90  % 


too  hieh 
cover 

too  low 
int.val 

too  hieh 
cover 

too  low 
int.val 


19(.06) 
262(.87) 

19(.061 
(.84,.91) 

11(.04) 
278  .93) 

11  .04) 
(.89.96) 


23(.08) 
243(.S1) 

34(.ll) 
(.77,.85j 

11(.04) 
278  .93) 

11  .04) 
(.90.93) 


27(.09) 

245(.S2) 
28(.09) 

(.77,.S6) 


14(.05) 
!74{.91) 

12(.04j 
.88..95) 


53(.17 

235  .78 
12(.04 

(.74,.83; 

32(.ll 

262(.87 
6(.02) 

(.84..9I 


46(.15 

233  .78 
2K.07 
(.73;.S2) 

23(.08) 
266(.89) 

IK. 04) 
(.85.92) 


2.0      80  % 


(.3857) 


too  hieh 
cover 

too  low 
int.val 

too  hieh 
cover 

too  low 
int.val 


18(.06) 
260(.87) 

22(.07j 
(.83;.91) 

9(.03) 
277(.92 

14(.05 
(.89.95 


31(.10) 
233  (.78 

36(.12 
(.73,.82) 

9(.03) 
269(.90] 

22  .071 
(.86..93) 


27(.091 
245  .82] 

28(.09J 
(.77;.86) 

14(.05] 
274C.91 

12(.04) 
(.88,-951 


47{.16) 
240  .80) 

13  .04) 
(.75,.85) 

22(.07) 
270(.90) 

8(.03) 
(.87..93) 


31 

243 
26 

(.77 


.10] 

.81 

.09^ 

.851 


15(.05) 

268(.89) 

17  .06) 

(.S6..93) 


3.0      80  % 


(.2399) 


90  % 


too  hieh 

cover 

too  low 

int.val 

too  hieh 

cover 

too  low 

int.val 


20(.07) 
262(.S7) 

18(.06) 
(.84,.9l) 

13(.04) 
283(.94) 
47.0 1J 

(.92.-97) 


20(.071 

247(.82 

33(.ll 

(.78,-87) 

13(.04) 
269  .90 

18f.06) 
(.86.93) 


27(.09) 

245(.82) 

28  .09) 

(.77,.86) 

14(.05) 
275  .92) 

11  .04) 
(-S9..95) 


36C12) 
250(.S3) 

14  .05) 
(.79,-88) 

20(.07 

274  .91 


6(.02) 
.88.-95) 


25(.08) 
240(.80 

35  .12 
(.75.-85) 

11(.04) 

267  .89) 

22  .07) 

(.85,.93) 


4.0      SO  % 


(.1492) 


too  hieh 
cover 

too  low 
int.val 

too  hieh 

cover 

too  low 

int.val 


14(.05) 
278(.93) 

87.03) 
(.90,.96) 

6(.02) 
286(.95) 

87.03) 
(.93,.98) 


14(.05) 
253  .84) 

33  .11) 
(.80;.S8) 

6{.02) 
26I(.87) 

33(.ll) 
(.83,.91) 


27(.09) 
245(.82) 

2S  .09) 
(.77;.86) 

14(.05) 

275(.92) 

11  (.04) 
(.S9..95) 


35  .12) 
24S(.S3) 

17(.06) 
(.78;.87) 

19(.06) 
274(.91) 

71-02) 
(.88,-95) 


20(.07 
239(.S0 

41  .14 
(.75;.84) 

10(.03) 
262  .87) 

28(.09) 
(.S4..91) 
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TABLE  2 

LENGTH  OF  C.I  FOR  In  P(D>  t}  (EXPONENTIAL  MODEL} 
N  =  50,  6  =  0.5,  p/=  1.0,  p2=  10.0 


time 

C.I  (%) 

BIN 

NOR 

m'le 

J.K 

B.T 

0.5 

80  % 

0.2161 

(0.033) 

0.1902 

(0.032) 

0.0889 

(0.012) 

0.1546 

(0.080) 

0.1267 

(0.038) 

90  % 

0.2709 
(0.042) 

0.2447 
(0.041) 

0.1141 
(0.016) 

0.2048 
(0.106) 

0.1619 
(0.047) 

1.0 

80  % 

0.3188 
(0.049) 

0.2882 

(0.042) 

0.1768 

(0.024) 

0.2022 
(0.066) 

0.1810 

(0.032) 

90  % 

0.4008 
(0.056) 

0.3717 
(0.055) 

0.2271 

(0.032) 

0.2680 

(0.088) 

0.2344 
(0.041) 

1.5 

SO  % 

0.4164 
(0.056) 

0.3798 
(0.053) 

0.2648 
(0.037) 

0.2822 
(0.088) 

0.2623 

(0.045) 

• 

90  % 

0.5237 
(0.071) 

0.4919 

(0.079) 

0.3400 
(0.048) 

0.3740 
(0.117) 

0.3398 
(0.061) 

2.0 

80  % 

0.5238 
(0.082) 

0.4802 
(0.063) 

0.3528 
(0.049) 

0.3773 
(0.120) 

0.3522 
(0.065) 

90  % 

0.6584 
(0.100) 

0.6253 

(0.076) 

0.4530 

(0.064) 

0.5001 

(0.159) 

0.4583 
(0.086) 

3.0 

80  % 

0.7496 
(0.131) 

0.6925 

(0.125) 

0.5287 

(0.074) 

0.5846 

(0.188) 

0.5449 
(0.107) 

90  % 

0.9404 
(0.164) 

0.9173 
(0.178) 

0.6790 
(0.095) 

0.7749 
(0.249) 

0.7059 
(0.140) 

4.0 

80  % 

1.0498 
(0.260) 

0.9973 
(0.311) 

0.7047 
(0.099) 

0.7988 
(0.261) 

0.73S0 
(0.148) 

90  % 

1.3136 
(0.323) 

1.4207 

(0.737) 

0.9049 

(0.127) 

1.0590 
(0.346) 

0.9589 
(0.194) 
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TABLE  3 

COVERAGE  RATIO  (EXP  MODEL) 
N  -  20,  8  =  0.5,  px=  1.0,  p2=  10.0 


time  C.I  (%) 


coverage 


BIN 


NOR 


MLE 


J.K 


B.T 


0.5      80  % 


(.7866) 


90 


/O 


too  high 
cover 

too  low 
int.val 

too  high 
cover 

too  low 
int.val 


14(.05) 

272(.91) 

14(.05) 

(.87,-94) 

4(.01) 
282C.94) 

14(.05) 
(.91  .97) 


40(.13 
228  .76 

32  .11) 
C-71,.81) 

40(.13) 
246(.82) 

14(.05) 
(.78-86) 


24(.08 
251  .84 

25  .08 
(.79,-88) 

14(.05j 

276  .92 

10  .03 

(.89.95 


106(.35) 

174(.58) 

20(.07) 

(.52.64) 


140(.47 
160(.53 
O(.'OO) 
(.48,-59) 


78(.26 
211( 
ll(.04) 

f. 65-76) 


26)         99(33} 
70)       20U.67) 
0?.00) 

(.62-72) 


1.0      80  % 


(.6203) 


90  % 


too  hieh 
cover 

too  low 
int.val 

too  hieh 

cover 

too  low 

int.val 


17(.06) 
250  .83) 

33(.U) 
(.79,.88) 

7(.02) 
27^(.92) 

18(.06) 
(.S8-95) 


40(.13) 

227(.76) 

33(.ll) 

(-71,.81) 

17(.06) 
2651.88) 

18(.06) 
(.S5-92) 


24(.08) 
250  .83) 

26  .09) 
(.79,-88) 

14(.05) 
277  .92) 

9(.03) 
(.89,-95) 


72(.24) 
210(.70) 

18(.06 
(.65,.75 

46(.15) 
245(.82) 

9(.03) 
(.77.-86) 


66(.22) 

224(.75) 
10(.03) 

(.70..80) 


259(. 


11) 
.86) 

8C03) 
(.82,.90) 


1.5      80  % 


(.4891) 


90  % 


too  hieh 
cover 

too  low 
int.val 

too  hiah 
cover 

too  low 
int.val 


15(.05) 
265(.88) 

20(.07) 
(.85,-92) 

15(.05) 

278  .93 
7(.02Y 

(.88.-93) 


26(.09) 
226(.75) 

48(.16) 
(.70,-80) 

15(.05) 

265(.88) 

20  (.07) 

(.85-92) 


24(,08) 
250(.83) 

26(.09) 
(.79,-88) 

14(.05 
277(.92 

9C03) 
(.89,-95) 


48(.16) 
247(.82) 

5(.02V 
(.78,-87) 

33(.ll) 
266  .89) 

H.00Y 
(.85,.94) 


40(.13) 

236  .79) 
24(.08) 

(.74..S3) 

19(.06) 

266(.S9) 
15(.05) 

(.85-94) 


2.0      80  % 


(.3857) 


90  % 


too  hiah 

cover 

too  low 

int.val 

too  hiah 

cover 

too  low 

int.val 


11(.04) 
267(.89) 

22.07) 
(.85,.93) 

ll(.04 

280  .93 
9(.03) 
.91,-96) 


24(.08) 
220(.73) 

56(.  19) 
(.68,-78) 

1H.04) 
267(.S9 

22(.07J 
(.85,.93v 


24(.08; 
250  (.83 

26  .09 
(.79,-88] 

14(.05) 
277(.92) 

9(.03j 
.89,-95) 


50(.17) 

245(.83) 
51-02) 
(.77,-86) 

28(.09) 
271  (.90) 

l?.00) 
(.87,-94) 


27(.09) 
240(.80) 

33  .11) 
(.75,-85) 

10(.03) 

267  .89) 

23  .08) 

(.85-93) 


3.0      80  % 


(.2399) 


90  % 


too  hieh 
cover 

too  low 
int.val 

too  hiah 

cover 

too  low 

int.val 


23(.08) 
271  .90) 

61-02) 
(.87-94) 

8(.03J 
286(.9^) 

6(.02) 
(.93-98) 


23(.08) 
249(.S3 

28(.09' 
(.79,-87 

8(.03) 
264(.88) 

28  .09) 

(.S4..92) 


23(.08] 


25K.84 
26(.09j 

(.79-S8; 

14(.05 

277  .92 
9(.03) 
(.89-95) 


40(.13) 
252  .84) 

8(.03) 
(.80,-88) 

22(.07) 
276  .92) 

2(.01) 
(.S9..9D) 


19(.06) 
238  .79) 

43(.14) 
(.75,-84) 

9(.03) 

260(.87 
31 

(.S3 


.10 
91) 


4.0      80  % 


(.1492) 


90  % 


too  hieh 
cover 

too  low 
int.val 

too  high 

cover 

too  low 

int.val 


19(.06) 
271  .90) 

10(.03) 

(.S7,.94) 

3(.01) 

287(.96) 
6(.02J 

(.93,-98) 


19(.06) 
224(.75) 

57(.19) 
(.70-80) 

3(.01) 
240(.8O) 

28(.09) 

(.75,.85) 


23(.08) 

36l 

.12) 

251  .84 

2531 

.84) 

26  .09) 

111 

'.04) 

(.79-88) 

(.81 

.89) 

14(.05 

277  .92 
9(.03) 
(.89-95) 


21(.07) 

277  .92) 
2(.01) 
(.89-95) 


15(.05 

231  .77 

54(.18 

(.72;.82) 

7(.03) 
26I(.87) 

32  .11) 
(.83.91) 
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TABLE  4 


LENGTH  OF  C 


C.I  FOR  In  P{D<t}  (EXPONEN1 
N  =  20,  6  =  0.5,  p.  =  1.0,  p2=  10.0 


XPONENTIAL  MODEL) 


1 

time 

C.I  (%) 

BIN 

NOR 

MLE 

J.K 

B.T 

0.5 

80  % 

0.3696 

(0.092) 

0.3026 

(0.092) 

0.1446 

(0.033) 

0.2083 

(0.079) 

0.1854 

(0.059) 

90  % 

0.4598 
(0.115) 

0.3897 
(0.122) 

0.1858 
(0.042) 

0.2713 
(0.103) 

0.2449 
(0.081) 

1.0 

80  % 

0.5451 

(0.131) 

0.4703 
(0.125) 

0.2877 
(0.066) 

0.3072 
(0.108) 

0.2899 
(0.081) 

90  % 

0.6797 
(0.164) 

0.6132 

(0.168) 

0.3695 
(0.085) 

0.4001 

(0.140) 

0.3816 
(0.110) 

1.5 

80  % 

0.7156 
(0.167) 

0.6320 
(0.159) 

0.4308 
(0.099) 

0.4668 
(0.180) 

0.4361 
(0.130) 

90  % 

0.8928 
(0.208) 

0.8351 

(0.227) 

0.5533 

(0.127) 

0.6079 
(0.235) 

0.5790 
(0.181) 

2.0 

SO  % 

0.9150 
(0.273) 

0.8328 
(0.306) 

0.5739 
(0.132) 

0.6452 
(0.258) 

0.5971 
(0.186) 

90  % 

1.1412 

(0.340) 

1.1847 
(0.837) 

0.7370 
(0.169) 

0.8403 
(0.336) 

0.7927 
(0.259) 

3.0 

80  % 

1.3762 
(0.567) 

1.3689 

(0.878) 

0.8601 

(0.198) 

1.0177 
(0.419) 

0.9324 

(0.307) 

90  % 

1.7105 
(0.663) 

2.2324 
(1.773) 

1.1045 
(0.255) 

1.3253 
(0.546) 

1.2427 
(0.429) 

4.0 

80  % 

2.1097 
(1.197) 

2.4656 
(2.060) 

1.1463 
(0.264) 

1.3966 

(0.584) 

1.2722 
(0.426) 

90  % 

2.5832 
(1.308) 

3.9914 

(2.720) 

1.4720 
(0.339) 

1.8188 
(0.761) 

1.6956 

(0.596) 
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TABLE  5 


CO\ 

N 


'ERAGE  RATIO  (EXP  MODEL) 
=  50,6=2/3,  pj  =  1.0,  P2=  10.0 


time  C.I  (%) 


coverage 


BIN 


NOR 


MLE 


J.K 


B.T 


0.5      80 


(.7231) 


90 


too  hiah 
cover 

too  low 
int.val 

too  hiah 
cover 

too  low 
int.val 


28(.09) 

249(.S3) 
23(.08) 

(.79,.S7) 

13(.04) 
277  (.92) 

10(.03) 
(.89;.95) 


50(.17) 

227(.76) 
23  (.08) 

20(.07] 

262(.S7^ 

10(.03 

(.84.91" 


27(.09; 
248  .83 

25  (.08 
(.78,.87j 

12(.04) 
273  .91 

15(.05; 
(.8S;.94 


173(.58) 
12H.40) 

61.02) 
(.35,.46) 

149(.50) 
148(.49) 

3(.01) 

(.44,.55) 


246(.82) 
54(.18) 
0(.00) 

(-14,.22) 

210(.70) 
90(.30) 
0(.00j 

(.25,.35) 


1.0      80  % 


(.5241) 


90  % 


too  high 
cover 

too  low 
int.val 

too  hieh 

cover 

too  low 

int.val 


20(.07 
251  (.83 


29(.10 

(.79..S8) 

13(.04) 
279(.93) 

8(.03V 
(.90..96) 


32(.ll) 

239  .80) 

29(.10) 

(.75,.84) 

20(.07) 

27U.90) 

9C03V 

(.S7..94) 


27(.09) 

248  .83 
25  .08< 
(.78,.87j 

11(.04) 

274  .91) 
15  (.05) 

(.S8..95) 


13K.44) 
160(.53) 

9(.03) 
(.48,.59) 

96(.32) 
200(.67) 

4(.01) 
(.61,-72)      ( 


150(.50) 

146(.49) 

47.01) 

(.43,-54) 

10H.34) 
198(.66) 

l(.0O) 
.61-71) 


1.5      80  % 


(.3799) 


90  % 


too  hiah 

cover 

too  low 

int.val 

too  hiah 

cover 

too  low 

int.val 


30(.10) 
239(.80) 

31(.10) 

(.75,.84) 

10(.03) 
279  .93) 

1T(.04) 
(.90,-96) 


30(.10 
239(.80 

3U.10) 
(.75,-84) 

15(.05) 
267  .89) 

18(.06) 
(.85.93) 


27(.09) 
248(.83) 

25  .08) 
(.78,.87) 

11(.04) 

274  .91) 

15  .05) 

(.88.95) 


86(.29) 
206(.69) 

8{.03V 
(.63,.74) 

53(.18 
244(.81 


3(.on 

,77,.86) 


S8(.29) 
199(.66) 

13(.04) 
(.61,-72) 

5H.17) 
243(.81) 

6C02) 
(.77,.85) 


2.0      80  % 


(.2753) 


90°, 


/o 


too  hiah 
cover 

too  low 
int.val 

too  hieh 
cover 

too  low 
int.val 


19(.06) 
254  (.85 

27  .09 
(.81,.89; 

9C.03) 
275(.92 

16  .05 
(.89,.95) 


19(.06) 

237(.79 

44(.15 

(.74,-84) 

9(.03) 
265(.88) 

26  .09) 
(.85.92) 


27(.09) 
248(.83) 

25  .08) 
(.78,-87) 

1H.04) 
274(.91) 

15  .05) 
(.88,.95) 


61(.20) 
228  .76) 

11  .04) 
(.71,.S1) 

34(.ll) 
261  .87) 

5(.02) 
(.83.-91) 


55(. 181 
227(.76 

18  .06 
(-71,-81) 

24(.0S; 
262(.S7< 

14(.05: 
(.84..91 


3.0       80  % 


(.1446) 


90  % 


too  hiah 
cover 

too  low 
int.val 

too  hiah 
cover 

too  low 
int.val 


32(.ll) 

252(.84) 

16(.05) 

(.80..88) 

16(.05 

277  .92 
7(.02) 
(.89,-95) 


32(.ll) 
220(.73) 

4S(.16) 
(.68.-78) 

6(.02) 
246(.82) 

48(.16) 
(.78.86) 


27(.09) 
248(.83 

25  .08 
(.78,.87 

12(.04) 

273  .91) 

15  .05) 

(.S8,.94) 


45.15) 

240  .80 

15  (.05 

(.75,.85 

19(.06) 

272  .91) 
91-03) 
(.87,-94) 


25(.08) 

239  (.80) 
36.12) 

(.75,.84) 

11(.04 

264  .88 
25  .08) 
(.84.92 


4.0      80  % 


(.0760) 


90  % 


too  hiah 
cover 

too  low 
int.val 

too  hiah 

cover 

too  low 

int.val 


35(.12) 
242  .81) 

23  (.08) 
(.76,.S5) 

12(.04) 
285(.95) 

3(.0l) 
(-93,-97) 


35(.12) 
196(.65) 

69(.23 
(.60,-71 

1(.00) 

276(.92) 

23(.08) 

(.89.95) 


27(.09; 

l         311 

.10) 

248  .83< 

\       25U 

'.84) 

25  .08' 

l         18i 

'.06) 

(.78.87) 

(.79 

.88) 

12(.04) 

273(.91) 

I5(.05j 

(.88.94) 


16(.05) 

274(.iJD 
10(.03 
(.88.95 


19(.06) 

233  (.78) 

4S(.16) 

(.73,'.82) 

8(.03) 
263  .88) 

29  .10) 
(.84,-91) 
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TABLE  6 

LENGTH  OF  C.I  FOR  In  P{D>t]  (EXPONENTIAL  MODEL) 
N=  50,  6  =  273,  P/-T.0,  P2=10.0 


time 

CI  (%) 

BIN 

NOR 

MLE 

J.K 

B.T 

0.5 

80  % 

0.2536 

(0.040) 

0.2261 

(0.038) 

0.1196 

(0.017) 

0.2232 
(0.092) 

0.1914 

(0.048) 

90  % 

0.3184 
(0.050) 

0.2913 
(0.049) 

0.1535 
(0.022) 

0.2958 
(0.122) 

0.2506 
(0.060) 

1.0 

80  % 

0.3848 
(0.056) 

0.3502 
(0.052) 

0.2383 
(0.034) 

0.2618 

(0.082) 

0.2510 

(0.048) 

90  % 

0.4840 
(0.070) 

0.4530 
(0.069) 

0.3060 

(0.044) 

0.3470 
(0.109) 

0.3243 
(0.064) 

1.5 

80  % 

0.5260 
(0.085) 

0.4823 
(0.080) 

0.3570 
(0.051) 

0.3841 
(0.134) 

0.3635 

(0.074) 

90  % 

0.6612 

(0.107) 

0.6282 

(0.107) 

0.4584 
(0.066) 

0.5091 

(0.178) 

0.4783 
(0.102) 

2.0 

80  % 

0.6815 

(0.126) 

0.6284 
(0.120) 

0.4757 
(0.068) 

0.5387 
(0.196) 

0.4977 
(0.106) 

90  % 

0.8556 

(0.157) 

0.8281 

(0.173) 

0.6108 

(0.088) 

0.7140 
(0.260) 

0.6584 
(0.144) 

3.0 

80  % 

1.0966 

(0.332) 

1.0598 

(0.468) 

0.7131 
(0.102) 

0.8739 

(0.332) 

0.7867 
(0.175) 

90  % 

1.3718 
(0.414) 

1.5369 
(0.930) 

0.9157 
(0.132) 

1.1583 
(0.440) 

1.0428 
(0.237) 

4.0 

80  % 

1.7374 
(0.773) 

1.9119 

(1.277) 

0.9505 
(0.137) 

1.2188 
(0.476) 

1.0833 
(0.248) 

90  % 

2.1575 
(0.901) 

3.1281 

(1.968) 

1.2205 
(0.176) 

1.6155 

(0.631) 

1.4384 
(0.336) 
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TABLE  7 

COVERAGE  RATIO  (EXP  MODEL) 
N  =  20,  9  =  2/3,  pL=  1.0,  p2=  10.0 


time  C.I  (%) 


BIN 


NOR 


MLE 


J.K 


B.T 


0.5      80  % 


(.7231) 


90  0/- 


0 


too  hiah 
cover 

too  low 
mt.val 

too  hiah 

cover 

too  low 

mt.val 


21(.07) 
257(.86) 

22(.07) 
(.82;.90) 

3{.01) 
282(.94) 

15(.05) 
(.91..97) 


50(.17) 
228  (.76) 

22(.07 
(-71,-81 

21(.07) 
264.88) 

15(.05) 
(.84.92) 


33(.ll) 

233(.78) 

34(.lli 

(.73,.82) 

1S(.06) 

268(.89) 
14(.05j 

(.S6..93) 


140(.47) 
155(.52) 
5(.02) 
(.46,.57) 

103(.34) 
194C.65) 


3(.01 
.59,7 


0) 


222(.74) 
77(.26) 
l(.OO) 
(-21,-31) 

187(.62) 

113(.38) 

0(.00) 

(.31.43) 


1.0      80  % 


(.5241) 


90  % 


too  hiah 
cover 

too  low 
int.val 

too  high 
cover 

too  low 
int.val 


31(.10) 
23K.77) 

38(.13) 
(.72,-82) 

14(.05) 
270(.9O) 

16(.05) 
(.87,-93) 


31(.10) 

231(77) 

3S(.13) 

(.70,.82) 

32(.ri) 

253(.84) 

16(.05) 
(.S0..S8) 


11) 

:^j(.78 
34U1 

(.73,.82) 

18(.06) 
268(.89) 

14  .05). 
(.S6..93) 


10U.34) 

193(.64) 

6(.02) 

(.59,-70) 

69(.23) 
230(.77) 

H.00) 
(.72,-81) 


112(.37) 

176(.59) 

12(.04) 

(.53..64) 

76(.25) 
217(,72) 

.  7(.02J 
(.67.-77) 


1.5      80  % 


(.3799) 


90  % 


too  hiah 

cover 

too  low 

int.val 

too  high 

cover 

too  low 

mt.val 


32(.ll] 

232(77 


36(.12) 
73;.82) 


14(.05 
273(.91 

13(.04 
(.S8..94) 


32(.ll) 

199(.66) 

69(.23) 

t.61,.72) 

14f.05) 
250(.83) 

36(-.12) 
(.79.881 


32U1) 

234(78) 

34(.ll) 

(.73;.83) 

18(.06 

268(.S9 
14(.05 
(.S6..93) 


8K.27) 
211(70) 

S(.03) 
(.65,76) 

52(.17) 
245(.82) 

3(.on 

(77..S6) 


67(.22) 
207(.69) 

261.09) 
(.64,-74) 

36(.12) 
24K.80) 

23(.08) 
(.76.85) 


2.0      SO  % 


(.2753) 


90  % 


too  hiah 
cover 

too  low 
int.val 

too  hiah 
cover 

too  low 
int.val 


21  (.07] 
25U.84 

28  .09 
(.79..S8; 

7(.02) 
287(.96) 

6(.02) 
(.93,.98) 


21  (.07) 
221(74) 

58(.19) 
(.69,79) 

7(.02) 
265(.S8) 

2S  .09) 
(.85.92) 


32(.ir 

235(781 

33(.ll 

(.74;.83) 

18(.06) 
269(.90) 

13(.04) 
(.86,-93) 


64(.21 
227  .76 

9(.03) 
(.71,-81) 

33(.ll 
262(.87 


(•8 


5(.'02) 
S3..9I) 


38(.13) 
227C76) 

35(.12 
(.71,-81 

15(.05) 
263(.88) 

22(.07) 
(.84.92) 


3.0      SO  % 


(.1446) 


90  % 


too  hiah 
cover 

too  low 
int.val 

too  hiah 
cover 

too  low 
int.val 


12(.04) 
274(.9l) 

14(.05) 
(.8S,.95) 

4(.01) 
282(.94) 

14(.05) 
(.91,-97) 


12(.04) 

226(.74) 

62  .21) 

(70,.80) 

4(.01 

234(7 
62(.21 

(.73..S3) 


32(.ll) 

33* 

(.74 


33(.ll) 
'  ..83) 


18(.06) 

270(.9() 
12  .04) 

{.87.. 93) 


50(.17) 
239  .80) 

11  .04) 
(.75).84) 

28(.09 

267(.S9 

5(.02) 

(.Si. 93) 


16(.05) 
237(79) 

47(.16) 
(.74;.84) 

9(.03) 
2637.88) 

2S(.09) 
(.84.-93) 


4.0      80  % 


(.0760) 


90  % 


too  hiah 
cover 

too  low 
int.val 

too  hiah 

cover 

too  low 

int.val 


21(.07) 
279(.93) 

0(.00) 
(.90,-96) 

4(.01) 

296(.98) 

Of  .00) 

(.97.-99) 


2H.07) 

215(,72) 

64(.21) 

(.67;.77) 

4(.01) 

232(77) 

64(.2lj 

(.73.82) 


18(.06) 
270(.90) 

12(.04) 
(■S7,.93) 


44(.15) 

244  .81) 

12(.04) 

(.77,.S6) 

24(.()8) 

27H.90) 

52-02) 

(.87,.94) 


15(.05) 
231  .77) 

54(.18) 
(-71,-82) 

6(.02) 
26IC-87) 

33(.ll) 
(.82.92) 


31 


TABLE  8 


LENGTH  OF  C 
N 


C.I  FOR  'nP(D>t}  (EXPONENTIAL  MODEL) 
N  =  20,  6  =  2,3,  p/=  1.0,  p2=  10.0 


time 

C.I  (%) 

BIN 

NOR 

MLE 

J.K 

B.T 

0.5 

80  % 

0.4304 
(0.109) 

0.3617 

(0.104) 

0.1981 

(0.053) 

0.3451 

(0.172) 

0.3168 

(0.128) 

90  % 

0.5359 
(0.173) 

0.4680 
(0.139) 

0.2544 
(0.068) 

0.4494 
(0.225) 

0.4302 
(0.183) 

1.0 

80  % 

0.6594 
(0.170) 

0.5788 
(0.161) 

0.3952 
(0.106) 

0.4630 
(0.207) 

0.4394 
(0.164) 

90  % 

0.8227 
(0.212) 

0.7618 
(0.225) 

0.5075 
(0.136) 

0.6030 

(0.298) 

0.5937 

(0.242) 

1.5 

80  % 

0.9426 

(0.323) 

0.8710 
(0.469) 

0.5931 
(0.161) 

0.7228 
(0.399) 

0.6560 

(0.256) 

90  % 

1.1756 

(0.402) 

1.2131 

(0.773) 

0.7616 

(0.207) 

0.9413 
(0.519) 

0.8865 

(0.379) 

2.0 

80  % 

1.2684 
(0.587) 

1.2540 
(0.899) 

0.8300 

(0.222) 

1.0229 
(0.603) 

0.9008 
(0.359) 

90  % 

1.5765 
(0.688) 

2.0258 
(1.817) 

0.9311 

(0.286) 

1.3322 
(0.785) 

1.2240 
(0.517) 

3.0 

80  % 

2.2188 
(1.311) 

2.5096 

(2.079) 

1.1951 
(0.400) 

1.6583 
(1.021) 

1.4301 

(0.574) 

90  % 

2.6995 
(1.400) 

4.1815 
(2.807) 

1.5347 
(0.914) 

2.1597 
(1.330) 

1.9457 
(0.830) 

4.0 

80  % 

3.5128 
(1.957) 

3.1967 

(2.724) 

1.6119 

(0.778) 

2.3080 
(1.443) 

1.9704 
(0.796) 

90  % 

4.0702 
(1.893) 

4.8194 
(3.183) 

2.0699 

(0.993) 

3.0057 
(1.880) 

2.6794 
(1.140) 

C.       ROBUSTNESS 

The  robustness  of  the  estimates  in  Chapter  II  was  studied  with  respect  to  an 
incorrect  model  assumption  about  the  distribution  of  the  sojourn  time  in  state  1.  In 
the  previous  simulations,  the  maximum  likelihood  estimator  used  the  known  correct 
model. 

The  data  for  the  simulation  experiment  in  this  section  is  generated  from  the 
following  three  state  semi-Markov  process:  Individuals  start  in  state  1  at  t  =  0.  The 
probability  of  a  transition  to  state  0  is  0,  and  to  state  2  is  1-6.  From  state  2,  the 
probability  of  a  transition  to  state  1  is  1.    State  0  is  an  absorbing  state.    The  sojourn 
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time  in  state  2  is  exponential  with  mean  l/p2.  The  sojourn  time  in  state  1  is  the  sum  of 
two  independent  exponentials  with  means  1/p,  and  1/p,;  that  is,  the  sojourn  time  in 
state  1  has  a  hypoexponential  distribution.  The  same  Fortran  program  for  the 
simulation  is  used  but  slightly  modified  for  the  above  change.  The  data  generated  are 
analyzed  by  the  same  Fortran  programs  for  each  estimator  as  in  the  first  section.  In 
particular,  the  maximum  likelihood  estimator  assumes  the  sojourn  time  in  state  1  has 
an  exponential  distribution  rather  than  the  true  hypoexponential  distribution. 

For  the  first  simulation  results  reported  in  Tables  9  to  12,  parameter  values  of  p, 
=  2.0,  p2  =  10.0,  p3  =  2.0,  and  8  =  0.5  are  used.  Again  two  different  numbers  of 
individuals  are  used;  they  are  20  and  50.  The  simulation  is  replicated  300  times,  the 
coverage  numbers  and  the  average  lengths  of  the  confidence  intervals  are  computed. 
The  actual  value  of  the  survival  function  is  computed  by  inverting  the  Laplace 
transform  of  the  passage  time  to  state  0  for  the  semi-Markov  process.  Some  computed 
values  of  the  survivor  function  can  be  found  in  Appendix  B. 

Table  9  reports  coverage  results  for  the  case  N  =  50  individuals.  The  use  of  the 
incorrect  model  for  the  maximum  likelihood  estimator  leads  to  the  majority  of  the 
intervals  being  too  low  for  t  =  0.5  and  t=  1.0.  The  binomial  confidence  interval  has  a 
slight  tendency  to  overcover.  The  confidence  interval  based  on  the  asymptotic  renewal 
estimator  undercover  for  t  =  0.5  and  t  =  1.0.  For  larger  values  of  t  they  have  the 
correct  coverage.  The  results  for  the  jackknife  and  bootstrap  individuals  are  similar  to 
those  in  Table  1.  Table  10  reports  the  average  lengths  of  the  confidence  intervals.  The 
average  lengths  for  the  maximum  likelihood  procedure  are  similar  to  those  of  Table  2. 
For  the  other  procedures,  the  average  lengths  are  smaller. 

Results  of  the  simulation  experiment  for  N  =  20  individuals  are  given  in  Table 
11  and  12.  Table  11  shows  the  coverage  results.  The  maximum  likelihood  procedure 
undercovers  for  t  =  0.5  and  t=  1.0  and  overcovers  for  t=  1.5  and  t  =  2.0.  The  binomial 
procedure  once  again  tends  to  overcover.  The  procedures  based  on  the  asymptotic 
renewal  estimator  have  the  correct  coverage  for  all  t  except  t  =  0.5.  The  average  lengths 
of  the  intervals  are  given  in  Table  12. 

Tables  13  to  16  report  results  of  simulation  of  a  semi-Markov  model  with  0=  2/3 
but  with  the  other  parameters  the  same;  Pj  =  2.0,  p2=  10.0,  p3  =  2.0  . 

Results  for  N  =  20  individuals  appear  in  Tables  15  and  16.  The  confidence 
intervals  based  on  the  asymptotic  renewal  estimator  have  the  correct  coverage  for  all  t 
except  0.5  and  1.0.  The  results  are  similar  to  chose  of  Tables  11  and  12. 
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TABLE  9 


COVER-AGE  RATIO  (HYPOEXP  MODEL) 
N  =  50,  6  =  0.5,  p2  =  2.0,  p2  -  10.0,  p3  =  2.0 


time  C.I  (%) 


length 


BIN 


NOR 


MLE 


J.K 


B.T 


0.5      80 


(.8652) 


90 


too  hieh 
cover 

too  low 
int.val 

too  high 
cover 

too  low 
int.val 


18(.06) 
26K.87) 

21(.0T 
(.83,-91 

10(.03) 
283(.94) 

7(.02) 
(.92,-97) 


35{.12) 
244  .81) 

21  .07) 
(.77;.86) 

18(.06 

275(.92 


7(.02j 
.89.-93 


0(.00) 
3(.0ij 

297(.99) 
(.00,-01) 

0(.00) 
6(.02) 

294(.98) 
(.00.-04) 


130(.43) 
165(.55) 
5(.02) 
(.49,.6I) 

99(.33) 
200(.67) 
1(.00) 
(.61,-72) 


149(.50) 
150(.50) 

(.44,-56) 

106(.35) 

194(.65) 

Ojf.OO) 

(.59..70) 


1.0      80  % 


(.6743) 


90  % 


too  hieh 
cover 

too  low 
int.val 

too  hieh 
cover 

too  low 
int.val 


22(.07) 
260(.87) 

18(.06) 
(.83,-91) 

11(.04 
280(.93 

9(.'03) 
(-9I..96) 


3(.01) 
229(.76 

33(.ll 
(.72,-81) 

22(.07 
269(.90 

9C03) 
(.S6..93) 


2(.01) 
145(.48) 

153(.51) 
(.43,-54) 

0(.00) 

2n(.7l) 

87(.29) 
(.66,.76) 


42(.14) 
242(.81 
16(.05< 

(.76,-85) 

28(.09) 
263  .88) 

9(.03V 
(.84.-91) 


39(.13) 

240  .SO) 
21  .07) 

(.75,-85) 

21(.07) 
265  .88) 

14(.05) 
(.85.-92) 


1.5      SO  % 


(.5157) 


90  % 


too  hish 
cover 

too  low 
int.val 

too  hieh 
cover 

too  low 
int.val 


:r(.o9) 

254(.S5) 

19(.06) 

(.81,-89) 

11(.04) 
276(.92) 

13(.04) 
(.89,-95) 


27(.09] 

235(.78 
38(.13 

(.74;.83; 

21(.07) 

261  .87 

18(.06 

(.83.91' 


6(.02) 

242(.81 
52(.17v 
(.76,-85) 

2761.93) 
20  .071 
(.S9..95 


3K.10) 
247  .82) 

22  .07 
(.78,-87) 

21(.07) 

271  .90) 
8(.02) 

(.87.-94) 


25(.08) 

243{.S1) 

32  .11) 

(.77;.85) 

19(.06) 
262  .87) 

19(.06) 
(.84..91) 


2.0      80  % 


(.3930) 


90  % 


too  hieh 
cover 

too  low 
int.val 

too  hiah 
cover 

too  low 
int.val 


25(.08] 
253  .84 

22  .07] 
(.80,-88] 

16(.05) 
274(.91) 

10(.03) 
(.88.-95) 


25(.08) 
240(.80) 

35.12) 
(.75;.85) 

16(.05] 

262(.87; 

22(.07) 

(.84.9T 


13(.04) 
264(.88) 

23(.08) 
(.84,.92) 

7(.02) 
282(.94) 

11(.04) 
(.91,-97) 


29(.10) 
250(.83) 

21  (.07) 
(.79,-88) 

20(.07) 
269  .90) 

11  .04) 
(.86,-93) 


24(.0S) 
239  .80) 

37.12) 
(.75,-84) 

12(.04) 
264  .88) 

24  .08) 
(.84.92) 


3.0      80  % 


(.2278) 


90  °, 


,  0 


too  hieh 
cover 

too  low 
int.val 

too  hieh 
cover 

too  low 
int.val 


25(.08] 
247(.82 

2S(.09 
(.78,-87) 

17(.06) 
273  .91) 

10  .03) 
(.88.94) 


25(.08) 
220(.73) 

55(.  18) 
(.68,-78) 

17(.06) 
255(.85 

2S(.()9 
(-81..89 


34(.1L 
253(.841 

13(.04] 
(.80,-88) 

20(.07 

275  .92 
5(.02j 
(.89.-93) 


31(.10) 
24S(.S3) 

21(.07) 
(.78,-87) 

18(.06) 
271(.90 

11  (.04 

(.87.94 


21( 

.07) 

241  ( 

.80) 

38( 

.13) 

(.76, 

.85) 

H( 

.04) 

266 

.89) 

23( 

.08) 

(.85-' 

.92) 

4.0      80  % 


(.1320) 


90  °/< 


too  hieh 
cover 

too  low 
int.val 

too  hieh 

cover 

too  low 

int.val 


11 
266 


.04 

:.89 

.08 
(.85,'.92) 

6(.02) 
285(.95) 

9(.03) 
(.93,-97) 


11(.04) 
235  .78 

54(.18 
(.74,-83) 

6(.02) 
271  .90) 
23(.08) 

(.87,.94) 


47(.16 
246  .82 

7(.02) 
(.78..S6) 

29(.10 

267(.S9 

4(.on 

(.85,-93) 


31  .10) 
244  .81) 

25(.08) 
(.77,.S6) 

19(.06) 
270(.90) 

11  (.04) 
(.87,.93) 


20(.07) 
23S  .79) 

42(.14) 
(.75,-84) 

10(.03) 

264(.S8) 

26(.09J 

(.84,.92) 
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TABLE  10 


LENGTH 

N 


H  OF  C.I  FOR  In  P(D>  tHHYPOEXP  MODEL) 
N  =  50,8  =  0.5,p1  =  !2.0,p2=10.0,p3  =  2.0 


time 

C.I  (%) 

BIN 

NOR 

MLE 

J.K 

B.T 

0.5 

80  % 

0.1688 

(0.029) 

0.1442 
(0.029) 

0.0886 
(0.01 1) 

0.0913 

(0.042) 

0.0806 

(0.022) 

90  % 

0.2110 
(0.037) 

0.1852 
(0.037) 

0.1138 
(0.014) 

0.1210 
(0.056) 

0.1043 
(0.028) 

1.0 

80  % 

0.2843 
(0.043) 

0.2553 
(0.041) 

0.1764 
(0.022) 

0.1666 
(0.053) 

0.1534 
(0.028) 

90  % 

0.3572 
(0.054) 

0.3291 
(0.053) 

0.2265 

(0.028) 

0.2208 
(0.071) 

0.1988 
(0.037) 

1.5 

80  % 

0.3958 
(0.058) 

0.3606 

(0.054) 

0.2641 

(0.033) 

0.2692 

(0.083) 

0.2496 
(0.047) 

90  % 

0.4978 
(0.073) 

0.4666 
(0.072) 

0.3392 
(0.042) 

0.3568 
(0.110) 

0.3240 
(0.064) 

2.0 

80  % 

0.5139 

(0.080) 

0.4710 
(0.075) 

0.3519 

(0.044) 

0.3777 
(0.117) 

0.3506 
(0.067) 

90  % 

0.6461 

(0.100) 

0.6130 
(0.100) 

0.4519 
(0.056) 

0.5006 
(0.155) 

0.4561 
(0.093) 

3.0 

80  % 

0.7840 
(0.160) 

0.7261 

(0.157) 

0.5274 
(0.066) 

0.5998 

(0.188) 

0.5576 
(0.111) 

90  % 

0.9832 
(0.199) 

0.9684 

(0.247) 

0.6773 
(0.085) 

0.7950 
(0.249) 

0.7261 
(0.153) 

4.0 

80  % 

1.1592 
(0.334) 

1.1276 
(0.480) 

0.7030 
(1.088) 

0.8242 
(0.260) 

0.7650 
(0.154) 

90  % 

1.4496 
(0.415) 

1.6667 
(1.017) 

0.9027 
(1.113) 

1.0924 

(0.345) 

0.9980 
(0.214) 
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TABLE  11 

COVERAGE  RATIO  (HYPOEXP  MODEL) 
N  =  20,6  =  0.5,  Pj  =  2.0,  p2*  10.0,  p3  =  2.0 


time  C.I  (%) 


length 


BIN 


NOR 


MLE 


J.K 


B.T 


0.5 


(.8652) 


/O 


90  % 


too  hiah 
cover 

too  low 
int.val 

too  hiah 
cover 

too  low 
int.val 


12(.04) 
278(.93" 
10(.03 

(.90,.96 

0(.00) 
290(.97) 

10(.03) 
(.95. .99) 


59(.20) 

231  .77) 
10  (.03) 

(.72..S2) 

37(.12 
260(.87 

3(.01) 
(.83,.92) 


0(.00) 
53(.18) 

247(.S2) 
(.13,.22) 

or.oo) 

109(.36) 
191(.64) 

(.31  .42) 


105(.35) 
183(.61) 


12(.04 

'.55, .67 


16 


)       ( 


73( 
220( 


.24 
.73 
7(.02) 
.68..7S) 


i35(.45) 
>2(.54) 
3(.01) 

48,-60) 


82(.37) 
217C.72) 

U.OOV 
(.67,-77) 


1.0      80  % 


(.6743) 


90  % 


too  hiah 
cover 

too  low 
int.val 

too  hi  ah 
cover 

too  low 
int.val 


16(.05) 
254(.85) 

30(.10) 
(.80,-88) 

5(.02) 
285(.95) 

10(.03l 
(.93.97) 


44(.15) 
226  .75) 

30  .10) 
(.70..80) 

16(.05) 
274C91) 

10  (.03) 
(.S8..95) 


1(.00) 
221(.74) 

78(.26) 
(.69,-79) 


0(.00) 

266(.S9) 

34(.ll) 

(.S5..92) 


45(.15) 
240(.80 
15f.05 

(.75,-85) 

29(.10) 
263  .88) 

8(.02)- 
(.S4..9I) 


32(.ll) 
24U.S0) 

27(.09) 
(.76,-85) 

11(.04 

271  .90 
18(.06 
(.87:.94) 


1.5      80  % 


(.5157) 


90  % 


too  high 
cover 

too  low 
int.val 

too  hieh 
cover 

too  iOW 
int.val 


2K.07) 
260(.S7) 

19(.06 
(.83,-91 

5(.02) 
276(.92) 

19(.06) 
(.89.95) 


21(.07) 

242(.S1) 

37.12) 

(.76,-85) 

21  (.07) 

260(.87) 

19(.06) 

/.83..91) 


S(.03) 
257(.86) 

35(.12) 

(.82,.90) 

21-01) 
283(.94 

15(.05 
(.92,.97 


40(.13) 
244(.81) 

16(.05 
(.77,.86' 

24(.08) 
273  .91) 

3(.01) 
(.S8..94) 


7(.02) 

264(.8S) 

29U0) 

(.84.92) 


2.0      80  % 


(.3930) 


90  % 


too  hiah 
cover 

too  low 
int.val 

too  hiah 
cover 

too  low 
int.val 


1U.04) 
270(.90) 

19(.06) 
(.87,-93) 

11  (.04) 
284(.95) 

5(.02) 
(-92,-97) 


28(.09) 
235(.7S 

37(.12^ 
(.74,-83; 

11  (.04 
270  .90 

19  .06) 
(.87.93) 


17(.06) 
260(.87) 

23(.08) 
(.83,-91) 

7(.02) 
284(.95) 


9(.03) 

,92,.97) 


39(.lbj 

246(.82; 

15(.05j 

(.78..861 

19(.06) 
277  .92) 

4(.01j 
(.89,.9o) 


18(.06) 

232  .77) 

50  .17) 

(.73,-82) 

7(.02) 

263(.S8) 
30  .10 
(.84..  91 


3.0       80  % 


(.2278) 


90  % 


too  hiah 
cover 

too  low 
int.val 

too  hiah 
cover 

too  low 
int.vai 


13  .04) 
273(.9l) 

14(.05) 
(.88,-94) 

5(.02) 
28U.94) 
14(.05) 

(.91. .96) 


13(.04) 
249  .83) 

3S  .13) 
(.79..87J 

5(.02) 
257(.86) 

38(.13) 
(.82.90) 


3K.10) 
254(.85) 

15  .05 
(.81,-89 

17(.06 
279(.93 

4(.01) 
(.90..96) 


37(.l2) 
251  .84) 

12  .04 
(.79,.88 

19(.06) 

277  .92) 

4(.01) 

(.S9..95) 


17C.06) 

230(.77 

53  .18 

(.72,-81 

6(.02) 
26i(.87 


j^.ll 
C.S3..91 


4.0       80  % 


(.1320) 


90  % 


too  hiah 
cover 

too  low 
int.val 

too  hiah 

cover 

too  low 

int.val 


4(.01) 
279(.93) 
17(.06) 

(.87,-94) 

4(.01) 
296C.98) 

0(.()()) 
(.97,.99) 


4  .01) 
226(.75 

70(.23 
(.70,-80 

2(.01) 

22S(.76) 

70(.23) 

(.71,-81) 


44(.15) 
245(.82) 

11  .04) 
(.79,-88) 

23(.08) 
273(.91) 

4C01) 
(.88,-94) 


38(.13) 
251  (.84) 

11  (.04) 
(.81,.89) 

19(.06) 

277(.92) 

4(:oij 

(.89,.95) 


17(.06) 

229  .76) 

54(.18) 

(.72,-82) 

6(.02) 
260(.87) 

34(.ll) 
(.S3;.91) 
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TABLE  12 


LENGTH  OF  C.I 


N 


3F  C.I  FOR  In  P(D>  t}(HYPOEXP  MODEL) 
=  20,6  =  0.5,pj  =  2.0,p2  =  10.0,p3  =  2.0 


time 

C.I  (%) 

BIN 

NOR 

MLE 

J.K 

B.T 

0.5 

80  % 

0.2891 

(0.077) 

0.2217 
(0.083) 

0.1438 
(0.030) 

0.1405 

(0.056) 

0.1296 

(0.046) 

90  % 

0.3591 
(0.095) 

0.2810 
(0.111) 

0.1846 
(0.038) 

0.1829 
(0.075) 

0.1750 
(0.066) 

1.0 

80  % 

0.4845 
(0.113) 

0.4131 

(0.107) 

0.2861 
(0.059) 

0.2694 
(0.105) 

0.2527 

(0.073) 

90  % 

0.6038 
(0.142) 

0.5365 

(0.143) 

0.3674 
(0.076) 

0.3509 
(0.136) 

0.3314 
(0.102) 

1.5 

80  % 

0.6810 

(0.175) 

0.5995 
(0.168) 

0.4284 
(0.088) 

0.4512 
(0.180) 

0.4190 
(0.128) 

90  % 

0.8496 

(0.218) 

0.7911 
(0.241) 

0.5501 
(0.113) 

0.5876 

(0.234) 

0.5517 

(0.174) 

2.0 

80% 

0.8843 
(0.230) 

0.7974 
(0.240) 

0.5707 
(0.118) 

0.6432 

(0.259) 

0.5944 
(0.290) 

90  % 

1.1030 

(0.286) 

1.0914 

(0.513) 

0.7329 
(0.151) 

0.8376 
(0.337) 

0.7864 
(0.256) 

3.0 

80  % 

1.4552 
(0.589) 

1.5826 
(1.296) 

0.8553 

(0.176) 

1.0345 
(0.419) 

0.9528 

(0.313) 

90  % 

1.8145 
(0.737) 

2.4881 
( 1.984) 

1.0983 

(0.227) 

1.3472 
(0.546) 

1.2605 
(0.427) 

4.0 

80  % 

2.3431 
(1.375) 

2.6289 
(2.141) 

1.1399 
(0.235) 

1.4289 
(0.581) 

1.3138 

(0.437) 

90  % 

2.8401 

(1.447) 

4.3272 
(2.778) 

1.4638 
(0.302) 

1.8610 

(0.757) 

1.7390 
(0.601) 
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TABLE  13 


OVERAGE  RATIO  (HYPOEXP  MODEL) 
N  =  50,6  =  2/3,  px  =  2.0,  p2=10.0.p3=2.0 


co\ 

N 


time  C.I  (%) 


leneth 


BIN 


NOR 


MLE 


J.K 


B.T 


0.5      80  % 


(.8214) 


90  % 


too  hieh 
cover- 
too  low 
int.val 

too  hieh 
cover 

too  low 
int.val 


25(.08) 
249  .83) 

26  .09) 
(.7S..S7) 

8(.03) 
284(.9<) 

8(.03) 
(.91.97) 


25(.08) 

249(.S3) 

26  .09) 

(.78..S7) 

25(.08 
267(.S9 

S(.03) 
(.8$,.93) 


0(.00) 

29S(.99) 
(.00..02) 


0(.00 
4(.01 

296(.9 
(.00..03 


1S6(.62) 

109(36) 

5p2J 

(.31,-42) 

157(.52) 
141(.47) 

2(.0D 
(.4I..53) 


220(.73) 
80(.27) 

0(.00V 
(•22,.32) 

194(.65) 
106(.35) 

Ojf.OO) 
(.30..41) 


1.0      80  % 


(.5788) 


90  % 


too  hieh 
cover 

too  low 
int.val 

too  hieh 
cover 

too  low 
int.val 


30(.10) 
254(.S5) 

16(.05) 
(.81,-89) 

10(.03 
2S0(.93, 

10  .03) 
(.91  .97) 


30(.10) 

244(.S1 
26(.09) 
(.77;.S6) 

15(.05) 
269  .90) 

16(.05) 
(.S6..93) 


0(.00) 
179(.60) 

121(.40) 
(.54.-65) 

or.oo) 

237(.79) 

63(.20 

(.74.84) 


65(.22) 
219(.73) 

1 1(.04) 
(.68,.78) 

39(.13) 

252(.84) 

9(.03V 

(.S0..SS) 


59(.20) 

225(.76) 

16(.05) 

(.7i;.si) 

32(.ll) 
261  .87 

7(.02V 
(.S3..90 


1.5      80  % 


(.3936) 


/O 


too  hieh 
cover 

too  low 
int.val 

too  hieh 
cover 

too  low 
int.val 


19(.06) 
253  .84) 

28  .09) 
(.80,-88) 

11(.04) 

274  .91) 

15  .05) 

(.S8..95) 


19(.06) 
242(.81) 

39  .13) 
(.76,-85) 

1H.04) 
261  .87) 

2S(.09) 
(.S3..91) 


S(.03) 
264L88) 

28(.09) 
(.84;.92) 

4(.on 

284(.95) 

12(.04) 
(.92.-97) 


42(.14) 
241  .80) 

17  .06} 
(.76,-85) 

17(.06) 

275  .92) 
8(.03j 
(.89..9D) 


29U0) 
240(.80) 

31(.10) 
(.75,-85) 

15(.05) 
265  .88) 

20(.07) 
(.85.92) 


2.0      80  % 


(.2652) 


/O 


too  hieh 
cover 

too  low 
int.val 

too  hieh 
cover 

too  low 
int.val 


26(.09; 
249  .83 

25  .08 
(.79,.87J 

5(.02) 
283(.94) 

12(.04) 
(.92,.97) 


26(.09: 
233  .78 


14(.05) 
261(.87) 

25(.08) 
(.S3,.91) 


27(.09) 
26U.87 

12(.04 
(.83,.91 

11(.04) 
285(.95) 

4ft)  1J 
(.93,-97) 


35(.12) 
249  .83 

16  .07 
(.79,-87 

1S(.06 

274(.91 


(.8 


S(.'03) 
8,-95 


24(.0S) 
235  .78) 

41  .14) 
(.74;.83) 

12(.04) 
262(.87) 

26  .09) 
(.84.91) 


3.0      80  % 


(.1197) 


90  % 


too  hieh 

cover 

too  low 

int.val 

too  hieh 
cover 

too  low 
int.val 


14(.05) 
271(.90) 

15(.05) 
(.87).94) 

5(.02) 

292(.97) 

3(.0O 

(.96,-99) 


14(.05) 
243  (.81) 

43  .14) 
(-77,-85) 

5{.02) 
252(.84) 
43(.14) 

(.80,-88) 


69(.23) 

227(.76) 

4(.01V 

(.7I,.8I) 

37(.12 
261  .87 

2(.0O 
(.83,-90 


35(.12) 

249  (.83) 

16  .05) 

(.79;.87) 

17(.06) 
275  .92) 

8(.03J 
(.S6,.9d) 


22(.07) 
229  .76) 

49  .16) 
(.72,-81) 

6(.02) 
26i(.S7) 

33  .11) 
(.83.90 


4.0      80  % 


(.0540) 


90  °/< 


too  high 
cover 

too  low 
int.val 

too  hieh 

cover 

too  low 

int.val 


13(.04) 
275(.92) 

12(.04) 
(.89,-95) 

4(.0O 
296L98) 

<){.!)()) 

(.97,-99) 


13(.04 
22K.74 

66(.22 
(.69,-79) 

4(.0D 
230(.77) 

66(.22) 
(-72.81) 


98(.33) 
199(.66) 

3(.0O 
(.61,-72) 

63(.2D 
237{.79) 

0(.00) 
(.74,-84) 


34(.ll 

248(.83 
IS  .06 
(.78,-87) 

15(.05) 

277C92) 


8(.03j 
.89,.9^) 


19(.06) 

229  .76) 

52(.17) 

(.72,-81) 

6(.02) 
260(.87) 

34.10 
(.83,-91) 
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TABLE  14 


LENGTH  OF  C.I  FOR  In  P[D>  t}(HYPOEXP  MODEL) 

N  =  50,e  =  2(3)p1  =  2.0,p2=10.0,P3  =  2.0 


time 

C.I  (%) 

BIN 

NOR 

MLE 

J.K 

B.T 

0.5 

80  % 

0.1956 

(0.032) 

0.1704 
(0.031) 

0.1194 

(0.014) 

0.1496 

(0.057) 

0.1351 

(0.035) 

90  % 

0.2450 
(0.040) 

0.2197 
(0.040) 

0.1533 
(0.0  IS) 

0.1983 
(0.076) 

0.1773 
(0.046) 

1.0 

80  % 

0.3463 

(0.049) 

0.3140 
(0.046) 

0.2378 
(0.028) 

0.2267 
(0.069) 

0.2140 
(0.042) 

90  % 

0.4355 
(0.062) 

0.4056 
(0.060) 

0.3054 
(0.036) 

0.3005 
(0.091) 

0.2768 
(0.053) 

1.5 

80  % 

0.5128 

(0.080) 

0.4700 

(0.075) 

0.3563 

(0.042) 

0.3751 
(0.012) 

0.3537 
(0.072) 

90  % 

0.6447 
(0.100) 

0.6117 

(0.100)  • 

0.4576 

(0.054). 

0.4972 
.  (0.159) 

0.45S8 
(0.095) 

2.0 

80  % 

0.7076 
(0.143) 

0.6531 

(0.135) 

0.4748 
(0.057) 

0.5405 

(0.174)- 

0.4977 
(0.106) 

90  % 

0.8881 

(0.176) 

0.8630 

(0.193) 

0.6097 
(0.073) 

0.7164 
(0.232) 

0.6584 

(0.144) 

3.0 

80  % 

1.2845 

(0.532) 

1.2267 
(0.526) 

0.7118 
(0.085) 

0.8833 

(0.287) 

0.8240 
(0.182) 

90  % 

1.5970 
(0.591) 

1.8844 
(1.190) 

0.9141 
(0.109) 

1.1708 
(0.380) 

1.0761 

(0.240) 

4.0 

80  % 

2.2724 
(1.102) 

2.6235 

(1.778) 

0.9488 
(0.113) 

1.2309 
(0.402) 

1.1452 
(0.259) 

90  % 

2.7889 
(1.216) 

4.1527 
(2.227) 

1.2184 
(0.145) 

1.6314 

(0.532) 

1.4950 

(0.336) 

39 


TABLE  15 

COVERAGE  RATIO  (HYPOEXP  MODEL) 
N  =  20,9  =  2/3,  Pj  =  2.0,  p2=10.0,p3  =  2.0 


time  C.I  (%) 


leneth 


BIN 


NOR 


MLE 


J.K 


B.T 


0.5      80  °A 


(.8214) 


too  hi2h 
cover 

too  low 
int.val 

too  hieh 
cover- 
too  low 
int.val 


277(.92) 

201.07) 
(.89,.95) 

31.01) 

2881.96) 

9(.03) 

(.94.98) 


26(.09] 
254C.S5 
20(.07j 

(.81,-88) 

26(.09 

265(.3S 

9(.03l 

(.S5..92) 


0(.00) 

54(.18) 

246(.82l 

(.14,-22) 

01.00) 
115(.38) 

1851.62) 
(.33,.44) 


124(.41) 
172(.57) 

41.01) 
(.52,.63) 

88(.29) 
211(.70) 

urn) 

(.65..76) 


198(.66) 
101(.34) 

11.00) 
(-28,-39) 

144(.48) 
1561.52) 
01.00) 
(.21.39) 


1.0      SO  % 


(.5788) 


90  % 


too  hieh 
cover 

too  low 
int.val 

too  hieh 

cover 

too  low 

int.val 


13(.04) 

267  .89) 

20(.07j 

(.S5,'.93) 


251.08) 
242(.81) 

33  .11 
(.76,.85 

22(.07) 

260  (.87 

181. 0.6 

(.S3..91 


0(.00) 

23S(.79 

62  .21 

(.75, .84 

0(.00) 
2741.91) 

26(.09j 
(.88,-95) 


62(.21) 
222(.74) 

161-05) 
(.69).79) 

37(.12) 
253(.84) 

10(.03) 
(.80..S8) 


57(.19) 
2161.72) 

27(.09) 
(.67;.77) 

20(.07) 
261  .87) 

19  .06) 
(.83.91) 


1.5      80  % 


(.3936) 


90  0/- 


0 


too  hieh 
cover 

too  low 
int.val 

too  hieh 

cover 

too  low 

int.val 


161.05) 
257(.86) 

27(.09) 
(.82..90) 

71.02) 
2831.94) 

10(.03) 
(.92.97) 


341.11) 
2181.73) 

481.16) 
(.68,.78) 

161.06) 
2571.86) 

271.09) 
(.S2.\90) 


111.04)      •  491.16) 


2641.88) 

25(.08) 
(.84,.92) 

41.01) 
2841.95) 

121.041 
1.92..97) 


2371.79) 
14(.05) 
(.74,-84) 

22(.07) 
269  .90) 

91.03) 
(.S6..93) 


251.08) 

23K.77) 

441.15) 

(.72,.S2) 

81.03) 
2641.88) 

2S  .09) 
(.84.92) 


2.0      80  % 


(.2652) 


90  % 


too  hieh 
cover 

too  low 
int.val 

too  hieh 
cover 

too  low 
int.val 


5(.02) 
2831.94) 

12(.04) 
(.92.97) 


161.05) 
2251.75) 

59(.20) 
(.70,-80) 

16(.05) 
245(.82) 

39(.13) 
(.77.-86) 


241.08) 
259(.86) 

17(.06) 
(.82,-90) 

tl(.04) 
2811.94) 

81. 03) 
(.9I..96) 


411-13) 
243(.8lj 

16  .07) 
(.77..85J 

201.07) 
2711.90) 

91.03V 
(.87..94) 


231.08; 
231  .77 

46  .15 
(.72;.82J 

71.02) 
2631. S8) 

30(.10) 
(.S3..92) 


3.0      SO  % 


(.1197) 


90  °A 


too  hi  ah 
cover 

too  low 
int.val 

too  hieh 
cover 

too  low 
int.val 


20(  .07) 
257(.86) 

23(.08 
(.82,-90 

31.01) 
2971.99) 

01.00) 
1.98..99) 


201.07) 
1891.63) 

91(.30) 
(.58,.68) 

31.01) 

2741.91 
23(.08 
1.S8..95) 


581.19) 
23U.77) 

11  .04) 
(.72,-82) 

26(.09 
269  (.90 

51.02) 
(.86,. 93) 


391. 13' 
2421.8 1 

19(.06) 
(.76,.S5) 

191.06) 

275  .92) 
61.02) 
(-S9..95) 


1S1.06) 
226(.75) 

56  .19) 
(.70,-80) 

61.02) 
2621.S7 

32  .11 
(.83,-92) 


4.0       80  % 


(.0540) 


90  % 


too  hieh 
cover 

too  low 
int.val 

too  hieh 

cover 

too  low 

int.val 


26(.09) 
274.91) 

0(  .00  V 
(.88,-95) 

61.02) 

294(.98) 

01.00J 

(.96,.99) 


61.02) 
196( 

98 
(.60, 

01.00) 
2021.67) 

98(.33) 
(.62.73) 


65(. 22 
226  .75 

91.03) 
(.71,-80) 

411.14) 
2561.85) 

31.01) 
(.8I..89) 


36(.12 
247  .82 

17  .06 
(.78,.87) 

16(.05 
278(.93 

(.90i.96) 


151.05 

229  .76 

56  .19 

(.72,-81) 

61.02) 
260(.87 

341.11 
(.S3,.91) 
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TABLE  16 

LENGTH  OF  C.I  FOR  In  P{D>  t}(HYPOEXP  MODEL) 
N  -  20,6  =  2/3,p1  =  2.0,p2=10.0,p3  =  2.0 


time 

C.I  (%) 

BIN 

NOR 

MLE 

J.K 

B.T 

0.5 

80  % 

0.3368 

(0.088) 

0.2706 

(0.088) 

0.1951 

(0.042) 

0.2539 
(0.149) 

0.2382 

(0.105) 

90  % 

0.4186 
(0.110) 

0.3466 
(0.119) 

0.2505 
(0.054) 

0.3306 
(0.194) 

0.3249 
(0.143) 

1.0 

80  % 

0.6008 

(0.158) 

0.5234 
(0.150) 

0.3892 

(0.084) 

0.3918 
(0.161) 

0.3675 
(0.124) 

90% 

0.7494 
(0.197) 

0.6858 

(0.208) 

0.4998 
(0.108) 

0.5103 

(0.209) 

0.4872 
(0.174) 

1.5 

80  % 

0.9084 
(0.310) 

0.8374 
(0.465) 

0.5840 

(0.129) 

0.6753 

(0.305) 

0.6142 

(0.213) 

90  % 

1.1331 

(0.387) 

1.1701 

(0.834) 

0.7500 
(0.166) 

0.8794 

(0.397) 

0.8144 

(0.298) 

2.0 

80  % 

1.3743 
(0.845) 

1.3323 
(1.078) 

0.7800 

(0.183) 

0.9871 

(0.470) 

0.8909 

(0.316) 

90  % 

1.6939 

(0.930) 

2.1678 

(2.013) 

1.0016 

(0.235) 

1.2855 
(0.612) 

1.1786 
(0.443) 

3.0 

80  % 

2.5888 
(1.497) 

2.9300 

(2.330) 

1.1776 

(0.357) 

1.6313 

(0.810) 

1.4596 

(0.534) 

90  % 

3.11S0 
(1.550) 

4J553 
(2.S63) 

1.5122 
(0.462) 

2.1245 
(1.055) 

1.9376 
(0.737) 

4.0 

80  % 

4.2251 
(2.031) 

3.3399 
(2.307) 

1.5894 
(0.763) 

2.2834 
(1.154) 

2.0356 
(0.749) 

90  % 

4.7925 
(1.883) 

4.5249 
(3.430) 

2.0410 

(0.980) 

2.9739 
(1.503) 

2.7125 
(1.03S) 
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IV.  CONCLUSIONS 

This  thesis  considers  the  problem  of  estimating  the  log  probability  for  a  semi- 
Markov  process  which  does  not  enter  a  particular  state  before  time  t. 

Simulation  is  used  to  study  procedures  for  obtaining  confidence  intervals  for  the 
In  probability  a  semi- Markov  process  enters  a  fixed  state  after  time  t.  The  data  arise 
from  observing  N  individuals.  Three  estimators  and  associated  confidence  interval 
procedures  are  compared:  The  three  estimators  use  different  amounts  of  information 
about  the  process.  The  maximum  likelihood  estimator  and  its  normal  confidence 
interval  uses  the  most  information.  An  estimator  based  on  the  observed  first  passage 
times  uses  the  least.  An  estimator  based  on  an  exponential  approximation  to  the 
survivor  function  of  the  first  passage  time  uses  an  intermediate  amount  of  information; 
confidence  intervals  for  this  last  estimator  are  obtained  using  jackknife  and  bootstrap 
procedures.    The  simulation  results  indicate  the  following. 

(1)  Larger  numbers  of  individuals  result  in  smaller  average  confidence  interval 
lengths. 

(2)  The  binomial  confidence  interval  tends  to  overcover.  This  is  true  in  general 
since  the  target  coverage  probability  is  a  lower  bound  on  the  true  coverage 
probability.  They  also  tend  to  have  larger  average  length. 

(3)  The  maximum  likelihood  confidence  intervals  have  the  correct  coverage  if  the 
model  on  which  they  are  based  is  correct.  If  the  model  is  incorrect  they  can 
either  overcover  or  undercover. 

(4)  If  a  jackknife  interval  does  not  cover  the  true  value  the  interval  will  tend  to 
be  'too  high'. 

(5)  The  confidence  intervals  using  the  asymptotic  renewal  estimator  have  the 
correct  coverage  for  largish  t. 

(6)  The  bootstrap  confidence  interval  requires  much  more  computation  than  the 
jackknife  confidence  interval.  Since  the  results  for  the  jackknife  and  bootstrap 
procedures  are  similiar,  ease  of  computation  appears  to  favor  the  jackknife 
procedure. 
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APPENDIX  A 
TRUE  PROBABILITY  TABLE  I 

TABLE  17 
EXPONENTIAL  MODEL  P(D>TIME}. 


time 

pL=  1.0  .p2=10.0 

Pl  =  1.0.p2  =  2.0 

Pj- 1.0^2*"  l0-° 

9  =  0.5 

0  =  0.5 

0  =  0.6667 

0.0 

1.0000 

1. 0000 

1.0000 

0.5 

0.7866 

0.7968 

0.7231 

1.0 

0.6203 

0.6503 

0.5241 

1.5 

0.4891 

0.5351 

0.3799 

2.0 

0.3S57 

0.4415 

*  0.2753 

2.5 

-  0.3042 

0.3646 

0.1995 

3.0 

0.2399 

0.3012 

0.1446 

3.5 

0.1S91 

0.2488 

0.1048 

4.0 

0.1492 

0.2055 

0.0~60 

4.5 

0.1176 

0.1698 

0.0551 

5.0 

0.0928 

0.1403 

0.0399 

5.5 

0.0731 

0.1159 

0.0289 

6.0 

0.0577 

0.0957 

0.0210 

6.5 

0.0455 

0.0791 

0.0152 

7.0 

0.0359 

0.0653 

0.0110 

7.5 

0.0283 

0.0540 

0.0080 

8.0 

0.0223 

0.0446 

0.005S 

8.5 

0.0176 

0.0368 

0.0042 

9.0 

0.0139 



0.0304 

0.0030 

The  sojourn  time  in  state  1  has  an  exponential  distribution  with  mean  l/p*.   The 
sojourn  time  in  state  2  has  an  exponential  distribution  with  mean  l/pr 
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APPENDIX  B 
TRUE  PROBABILITY  TABLE  II 

TABLE  18 
HYPOEXPONENTIAL  MODEL  P{D>TIME}. 


time 

p1  =  2.0  ,p2=10.0 

Pl  =  2.0,p2=2.0 

p1  =  2.0,p2=10.0 

p^  =  2.0  .0  =  0.5 

pr  =  2.0  .6  =  0.5 

p,  =  2.0.8  =  0.666'' 

0.0 

1.0000 

1.0000 

1.0000 

0.5 

0.8652 

0.8670 

0.8214 

1.0 

0.6743 

0.6897 

0.5788 

1.5 

0.5157 

0.5519 

0.3936 

2.0 

0.3930 

0.4465  * 

0.2652 

2.5 

0.2992 

0.3268 

0.1783 

3.0 

0.2278 

0.2952 

0.1197 

3.5 

0.1734 

0.2402 

0.0804 

4.0 

0.1320 

0.1954 

0.0540 

4.5 

0.1005 

0.1590 

0.0363 

5.0 

0.0765 

0.1293 

0.0244 

5.5 

0.0582 

0.1052 

0.0164 

6.0 

0.0443 

0.0856 

0.0110 

6.5 

0.0337 

0.0697 

0.0074 

7.0 

0.0257 

0.0567 

0.0050 

7.5 

0.0196 

0.0461 

0.0033 

8.0 

0.0149 

0.0375 

0.0022 

8.5 

0.0113 

0.0305 

0.0015 

9.0 

0.0086 

0.0248 

0.0010 

The  sojourn  time  in  state  1  is  the  sum  of  exponential  distribution  with  mean  1/p, 
and  mean  l/p3<  The  sojourn  time  in  state  2  has  an  exponential  distribution  with  mean 
1/P2. 


44 


APPENDIX  C 
AVERAGE  LENGTH  OF  C.I  FOR  P{D>T} 

Recorded  in  the  Appendix  are  the  average  length  of  confidence  intervals  for 
P{D>t}  corresponding  to  those  given  in  the  thesis  for  In  P(D>t}.    If  a  confidence  limit 
is  greater  than  1,  it  is  set  equal  to  1.    If  a  confidence  limit  is  less  than  0,  it  is  set  equal 
to  0.   The  results  are  similar  to  those  in  the  previous  tables. 


TABLE  IS 

LENGTH  OF  C.I  FOR  P[D>t}  (EXPONENTIAL  MODEL) 
N  =  50,  6  =  0.5,  p\  =  1.0,  p2=  10.0 


time 

C.I'(%) 

BIN 

NOR 

VILE 

J.K 

B.T 

0.5 

80  % 

0.1651 
(0.014) 

0.1469 
(0.014) 

0.0692 

(0.007) 

0.1259 
(0.058) 

0.1044 
(0.027) 

90  % 

0.2054 
(0.018) 

0.1887 
(0.019) 

0.0889 
(0.009) 

0.1669 

(0.077) 

0.1330 
(0.033) 

1.0 

80  % 

0.1915 
(0.006) 

0.1741 
(0.007) 

0.1078 
(0.007) 

0.1278 
(0.036) 

0.1154 

(0.015) 

90  % 

0.2389 
(0.008) 

0.2236 
(0.009) 

0.1384 
(0.009) 

0.1694 

(0.048) 

0.1483 
(0.018) 

1.5 

80  % 

0.1967 

(0.002) 

0.1795 
(0.003) 

0.1264 

(0.005) 

0.1381 

(0.035) 

0.1289 
(0.014) 

90  % 

0.2455 
(0.003) 

0.2305 

(0.003) 

0.1623 
(0.006) 

0.1830 
(0.047) 

0.1650 

(0.019) 

2.0 

80  % 

0.1913 

(0.007) 

0.1740 
(0.007) 

0.1320 

(0.002) 

0.1427 
(0.035) 

0.1334 
(0.015) 

90  % 

0.2387 
(0.008) 

0.2234 
(0.009) 

0.1695 
(0.002) 

0.1892 
(0.046) 

0.1711 

(0.019) 

3.0 

80  % 

0.1706 
(0.012) 

0.1526 

(0.013) 

0.1220 
(0.008) 

0.1325 
(0.031) 

0.1233 
(0.016) 

90  % 

0.2123 
(0.016) 

0.1960 
(0.016) 

0.1567 
(0.010) 

0.1756 
(0.041) 

0.1569 
(0.021) 

4.0 

80  % 

0.1459 

(0.017) 

0.1271 

(0.017) 

0.1007 

(0.013) 

0.1092 

(0.02S) 

0.1009 

(0.018) 

90  % 

0.1820 
(0.021) 

0.1631 

(0.023) 

0.1293 

(0.017) 

0.1447 
(0.037) 

0.1283 

(0.023) 
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TABLE  19 

LENGTH  OF  C.I  FOR  P[D>t}  (EXPONENTIAL  MODEL) 
N  =  20.  8  =  0.5,  Pi  =  1.0,  p:=  10.0 


time 

C.I  (%) 

BIN 

NOR 

MLE 

J.K 

B.T 

0.5 

80  % 

0.2689 

(0.037) 

0.2262 

(0.045) 

0.1111 

(0.01S) 

0.1746 

(0.065) 

0.1584 

(0.055) 

90  % 

0.3281 
(0.046) 

0.2891 
(0.060) 

0.1426 
(0.023) 

0.2274 
(0.085) 

0.2094 
(0.078) 

1.0 

SO  % 

0.3087 
(0.020) 

0.2702 
(0.023) 

0.1708 
(0.018) 

0.1913 

(0.050) 

0.1820 

(0.035) 

90  % 

0.3775 
(0.024) 

0.3469 
(0.030) 

0.2194 

(0.023) 

0.2492 
(0.065) 

0.2356 

(0.047) 

1.5 

80  % 

0.3169 
(0.010) 

0.2792 
(0.011) 

0.1985 
(0.010) 

0.2165 
(0.051) 

0.2030 
(0.029) 

90  % 

0.3876 
(0.013) 

0.3585 
(0.014) 

0.2549 

(0.013) 

0.2820 

(0.-067) 

0.2617 

(0.040) 

2.0 

80  % 

0.3082 
(0.019) 

0.2698 
(0.021) 

0.2058 
(0.007) 

0.2244 
(0.049) 

0.2077 
(0.028) 

90  % 

0.3768 

(0.024) 

0.3463 
(0.027) 

0.2643 
(0.009) 

0.2923 
(0.064) 

0.2650 

(0.037) 

3.0 

80  % 

0.2758 
(0.032) 

0.2345 
(0.037) 

0.1886 
(0.022) 

0.2027 
(0.049) 

0.1854 

(0.034) 

90  % 

0.3365 
(0.040) 

0.3001 

(0.049) 

0.2422 
(0.029) 

0.2640 
(0.063) 

0.2338 
(0.044) 

4.0 

80  % 

0.2377 
(0.045) 

0.1895 
(0.057) 

0.1554 

(0.032) 

0.1627 
(0.050) 

0.1485 
(0.040) 

90  % 

0.2896 

(0.055) 

0.2389 

(0.077) 

0.1992 

(0.042) 

0.2118 

(0.065) 

0.1860 

(0.049) 
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TABLE  20 

LENGTH  OF  C.I  FOR  P(D>t}  (EXPONENTIAL  MODEL) 
N  =  50,8-2/3,  Pj- 1.0,  p2=10.0 


time 

C.I  (%) 

BIN 

NOR 

MLE 

J.K 

B.T 

0.5 

80  % 

0.1719 
(0.012) 

0.1601 
(0.012) 

0.0855 

(0.008) 

0.1873 

(0.074) 

0.1645 

(0.044) 

90  % 

0.2160 
(0.015) 

0.2056 
(0.016) 

0.1098 
(0.010) 

0.2482 
(0.098) 

0.2158 
(0.056) 

1.0 

80  % 

0.1962 

(0.003) 

0.1790 
(0.004) 

0.1223 
(0.006) 

0.1523 

(0.042) 

0.1479 
(0.023) 

90  % 

0.2448 
(0.004) 

0.2298 
(0.005) 

0.1571 

(0.007) 

0.2019 

(0.056) 

0.1897 
(0.031) 

1.5 

80  % 

0.1911 

(0.007) 

0.1737 
(0.007) 

0.1320 
(0.002) 

0.1542 
(0.043) 

0.1474 
(0.019) 

90  % 

0.2383 
(0.009) 

0.2231 
(0.009) 

0.1694 

(0.002) 

0.2044 
(0.057) 

0.1907 
(0.025) 

2.0 

80  % 

0.1772 
(0.011) 

0.1595 
(0.012) 

0.1269 
(0.006) 

0.1496 

(0.042) 

0.1397 
(0.019) 

90  % 

0.2208 
(0.015) 

0.2048 
(0.016) 

.  0.1629 

(0.008) 

0.1982 

(0.055) 

0.1803 

(0.024) 

3.0 

80  % 

0.1435 
(0.019) 

0.1245 
(0.020) 

0.0996 

(0.013) 

0.1178 
(0.037) 

0.1073 
(0.022) 

90  % 

0.1780 
(0.024) 

0.1598 
(0.026) 

0.1279 
(0.017) 

0.1561 

(0.049) 

0.1372 
(0.028) 

4.0 

80  % 

0.1140 
(0.021) 

0.0932 

(0.024) 

0.0701 
(0.015) 

0.0812 
(0.031) 

0.0733 
(0.021) 

90  % 

0.1409 
(0.026) 

0.1185 
(0.032) 

0.0900 
(0.020) 

0.1076 

(0.041) 

0.0934 

(0.02S) 
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TABLE  21 

LENGTH  OF  C.I  FOR  P{D>t}  (EXPONENTIAL  MODEL) 
N  =  20,  6  =  2/3,  Pj-1.0,  p2=  10.0 


time 

C.I  (%) 

BIN 

NOR 

MLE 

J.K 

B.T 

0.5 

80  % 

0.2874 
(0.030) 

0.2472 
(0.033) 

0.1374 
(0.021) 

0.3199 
(0.230) 

0.2984 
(0.158) 

90  % 

0.3508 
(0.037) 

0.3167 
(0.045) 

0.1765 
(0.027) 

0.4166 
(0.300) 

0.4154 
(0.238) 

1.0 

80  % 

0.3151 
(0.012) 

0.2771 
(0.013) 

0.1926 

(0.014) 

0.2595 
(0.098) 

0.2520 

(0.077) 

90  % 

0.3854 
(0.015) 

0.3560 
(0.016) 

0.2473 
(0.018) 

0.3380 
(0.127) 

0.3355 
(0.121) 

1.5 

SO  % 

0.3057 
(0.022) 

0.2670 

(0.024) 

0.2046 
(0.011) 

0.2539 
(0.076) 

0.2387 
(0.056) 

• 

90  %- 

0.3736 
(0.027) 

0.3428 
(0.031) 

0.2627 
(0,014) 

0.3306 

(0.099) 

0.3084 
(0.078) 

2.0 

80  % 

0.2S47 
(0.034) 

0.2440 
(0.038) 

0.2947 
(0.021) 

0.2338 
(0.071) 

0.2140 
(0.051) 

90  % 

0.3475 
(0.042) 

.0.3123 
(0.051) 

0.2499 
(0.028) 

0.3044 
(0.093) 

0.2730 
(0.066) 

3.0 

80  % 

0.2334 
(0.047) 

0.1842 
(0.061) 

0.1515 
;0.036) 

0.1700 
(0.066) 

0.1542 
(0.054) 

90  % 

0.2842 
(0.056) 

0.2316 
(0.081) 

0.1940 
(0.048) 

0.2213 
(0.086) 

0.1936 
(0.066) 

4.0 

80  % 

0.1893 
(0.051) 

0.1253 
(0.075) 

0.1069 
(0.039) 

0.1117 
(0.057) 

0.1024 

(0.049) 

90  % 

0.2316 
(0.060) 

0.1536 
(0.096) 

0.1351 

(0.052) 

0.1455 
(0.074) 

0.12S1 
(0.060) 
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TABLE  22 

LENGTH  OF  C.I  FOR  PfD>t)  (HYPOEXP  MODEL) 
N  =  50,0  =  0.5,  pj  =  2.0,  p2=  10.0,  p3  =  2.0 


time 

C.I  (%) 

BIN 

NOR 

MLE 

J.K 

B.T 

0.5 

80  % 

0.1417 
(0.018) 

0.1217 
(0.019) 

0.0692 
(0.006) 

0.0813 
(0.035) 

0.0722 

(0.018) 

90  % 

0.1758 
(0.022) 

0.1574 
(0.024) 

0.0888 
(0.008) 

0.1077 
(0.047) 

0.0931 
(0.022) 

1.0 

80  % 

0.1855 

(0.009) 

0.1680 
(0.010) 

0.1078 

(0.007) 

0.1110 

(0.032) 

0.1024 

(0.014) 

90  % 

0.2313 
(0.012) 

0.2157 
(0.013) 

0.1384 
(0.009) 

0.1472 
(0.042) 

0.1317 
(0.019) 

1.5 

80  % 

0.1964 

(0.003) 

0.1792 
(0.003) 

0.1265 

(0.004) 

0.1344 
(0.035) 

0.1245 
(0.016) 

90  % 

0.2451 
(0.004) 

0.2301 

(0.004) 

0.1625 
.  (0.006) 

0.1182 

(0.047) 

0.1597 
(0.021). 

2.0 

80  % 

0.1920 
(0.006) 

0.1747 
(0.006) 

0.1322 

(0.002) 

0.1413 
(0.036) 

0.1308 
(0.017) 

90  % 

0.2395 

(0.008) 

0.2243 
(0.008) 

0.1698 

(0.002) 

0.18-72 
(0.048) 

0.1676 
(0.021) 

3.0 

80  % 

0.1675 
(0.014) 

0.1495 
(0.015) 

0.1223 

(0.007) 

0.1265 

(0.032) 

0.1174 
(0.019) 

90  % 

0.2185 
(0.018) 

0.1919 
(0.015) 

0.1573 
(0.009) 

0.1677 

(0.043) 

0.1496 
(0.022) 

4.0 

80  % 

0.1390 
(0.018) 

0.1198 
(0.019) 

0.1009 
(0.011) 

0.0988 
(0.028) 

0.0919 

(0.019) 

90  % 

0.1723 
(0.023) 

0.1537 
(0.025) 

0.1296 

(0.015) 

0.1310 

(0.037) 

0.1168 
(0.023) 
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TABLE  23 

RP(D>t;  i 

N  =  20,  9  =  0.5,  p  j  =  2.0,  p2  =  10.0,  p5  =  2.0 


LENGTH  OF  CI  FOR  P{D>  t}  (HYPOEXP  MODEL) 
=  20,  8: 


time 

C.I  (%) 

BIN 

NOR 

MLE 

J.K 

B.T 

0.5 

80  % 

0.2328 
(0.044) 

0.1839 

(0.057) 

0.1108 
(0.016) 

0.1260 

(0.051) 

0.1177 

(0.044) 

90  % 

0.2836 
(0.053) 

0.2314 
(0.077) 

0.1423 
(0.021) 

0.1641 
(0.067) 

0.1592 
(0.065) 

1.0 

80  % 

0.3005 

(0.023) 

0.2614 

(0.026) 

0.1710 
(0.016) 

0.1736 
(0.050) 

0.1633 

(0.032) 

90  % 

0.3671 
(0.029) 

0.3354 
(0.034) 

0.2196 
(0.021) 

0.2261 
(0.066) 

0.2105 

(0.043) 

1.5 

80  % 

0.3163 
(0.011) 

0.2785 
(0.011) 

0.1990 
(0.010) 

0.2110 

(0.053) 

0.1954 
(0.031) 

90  % 

0.3869 
(0.013) 

0.3576 
(0.015) 

0.2556 

(0.012) 

0.2747 
(0.069). 

0.2502 

(0.040) 

2.0 

80  % 

0.3104 
(0.016) 

0.2722 
(0.017) 

0.2067 
(0.006) 

0.2195 

(0.050) 

0.2017 

(0.032) 

90  % 

0.3796 
(0.020) 

0.3495 
(0.023) 

0.2654 

(0.007) 

0.2859 
(0.065) 

0.2571 
(0.038) 

3.0 

80  % 

0.2709 

(0.035) 

0.2289 
(0.040) 

0.1896 
(0.020) 

0.1917 

(0.047) 

0.1756 
(0.037) 

90  % 

0.3304 
(0.043) 

0.2924 
(0.054) 

0.2434 
(0.025) 

0.2497 
(0.062) 

0.2211 
(0.045) 

4.0 

80  % 

0.2276 
(0.046) 

0.1771 
(0.061) 

0.1560 
(0.029) 

0.1470 
(0.047) 

0.1354 
(0.040) 

90  % 

0.2774 
(0.055) 

0.2223 

(0.082) 

0.2001 

(0.038) 

0.1914 

(0.061) 

0.1696 
(0.050) 
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TABLE  24 

:tTH  OF  C.I  FOR  P{D>t}  (HYPOEXP  MO 
N  =  50,6  =  2,3,  p,  =  2.0,  P2=10.0,  p3  =  2.0 


time 

C.I  (%) 

BIN 

NOR 

NILE 

J.K 

B.T 

0.5 

80  % 

0.1560 

(0.015) 

0.1375 

10.016) 

0.0855 

(0.007) 

0.1349 

(0.050) 

0.1237 
(0.034) 

90  % 

0.1938 
(0.019) 

0.1766 
(0.020) 

0.1098 
(0.009) 

0.1788 
(0.066) 

0.1624 
(0.045) 

1.0 

80  % 

0.1943 
(0.005) 

0.1771 
(0.005) 

0.1226 

(0.005) 

0.1330 
(0.036) 

0.1257 
(0.017) 

90  % 

0.2425 
(0.006) 

0.2274 
(0.006) 

0.1574 
(0.006) 

0.1763 

(0.048) 

0.1615 
(0.021) 

1.5 

80  % 

0.1921 
(0.006) 

0.1748 
(0.007) 

0.1324 
(0.001) 

0.1435 

(0.038) 

0.1346 
(0.014) 

90  % 

0.2397 
(0.008) 

0.2245 
(0.00.8) 

0.1700 
(0.002) 

0.1901 
(0.050) 

0.1722 
(0.019) 

2.0 

80  % 

0.1747 
(0.013) 

0.1569 
(0.014) 

0.1273 
(0.005) 

0.1351 

(0.035) 

0.1257 
(0.015) 

90  % 

0.2176 
(0.017) 

0.2015 

(0.018) 

0.1634 

(0.006) 

0.1791 

(0.046) 

0.1606 
(0.020) 

3.0 

80  % 

0.1327 
(0.019) 

0.1132 
(0.021) 

0.0998 
(0.011) 

0.0957 
(0.028) 

0.0885 
(0.018) 

90  % 

0.1644 
(0.024) 

0.1451 
(0.028) 

0.1281 
(0.014) 

0.1268 
(0.037) 

0.1123 
(0.022) 

4.0 

80  % 

0.0997 
(0.022) 

0.0769 
(0.026) 

0.0700 
(0.013) 

0.0558 
(0.022) 

0.0544 
(0.016) 

90  % 

0.1230 
(0.027) 

0.0963 
(0.036) 

0.0898 
(0.016) 

0.0779 

(0.029) 

0.0690 

(0.021) 
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TABLE  25 

LENGTH  OF  C.I  FOR  P(D>t)  (HYPOEXP  MODEL) 
N  =  20,  0=2,3.  p,  =  2.0,  p2=  10.0,  p3  =  2.0 


time 

C.I  (%) 

BIN 

NOR 

MLE 

J.K 

B.T 

0.5 

80  % 

0.2554 

(0.040) 

0.2112 

(0.047) 

0.1370 
(0.018) 

0.2449 
(0.197) 

0.2360 

(0.138) 

90  % 

0.3113 
(0.049) 

0.2683 
(0.065) 

0.1759 
(0.023) 

0.3189 

(0.257) 

0.3274 
(0.192) 

1.0 

80  % 

0.3126 
(0.015) 

0.2746 
(0.016) 

0.1935 
(0.011) 

0.2194 

(0.076) 

0.2105 

(0.056) 

90  % 

0.3823 
(0.019) 

0.3525 
(0.021) 

0.2484 
(0.015) 

0.2857 
(0.098) 

0.2744 
(0.080) 

1.5 

80  % 

0.3084 
(0.020) 

0.2699 
(0.022) 

0.2066 
(0.007) 

0.2298 
(0.060) 

0.2134 
(0.040) 

90  % 

0.3770 

(0.025) 

0.3465 
(0.030) 

0.2653 
(0.009) 

0.2993 

(0.078) 

0.2716 

(0.050) 

2.0 

80% 

0.2799 
(0.039) 

0.2381 
(0.047) 

0.1971 
(0.017) 

0.2095 
(0.057) 

0.1915 

(0.042) 

90  % 

0.3417 
(0.048) 

0.3044 
(0.062) 

0.2530 
(0.022) 

0.2728 
(0.075) 

0.2403 
(0.051) 

3.0 

80  % 

0.2180 

(0.047) 

0.1648 
(0.064) 

0.1533 

(0.030) 

0.1414 

(0.052) 

0.1297 

(0.043) 

90  % 

0.2657 
(0.057) 

0.2054 
(0.085) 

0.1965 
(0.040) 

0.1842 
(0.068) 

0.1619 
(0.054) 

4.0 

80  % 

0.1692 
(0.048) 

0.0964 
(0.074) 

0.1076 
(0.033) 

0.0849 
(0.042) 

0.0794 
(0.037) 

90  % 

0.2080 
(0.056) 

0.1163 

(0.092) 

0.1361 

(0.045) 

0.1105 
(0.055) 

0.0995 
(0.046) 

52 


LIST  OF  REFERENCES 


1.  Gallagher,  R.  M.,  A  Simulation  Study  of  Esimate  of  a  First  Passage  Time 
Distribution  for  a  Censored  Semi- Markov  Process,  Masters  Thesis,  Naval 
Postgraduate  School,  Monterey,  California,  September  1986. 

2.  Park,  D.  K.,  Parametric  and  Nonparametric  Estimation  of  the  Mean  Numbers  in 
Service  for  an  MlGj  co  Oueue.  Masters  Thesis,  Naval  Postgraduate  School, 
Monterey,  California,  March  1986. 

3.  Jacobs,  P.  A..  Results  of  a  Simulation  Studv  of  Esimates  of  a  First  Passage  Time 
Distribution  for  Censored  Semi-Markov  Process,  Naval  Postgraduate  School, 
Monterey,  California,  19S6. 

4.  Bickel,  P.  J.  and  Doksum,  K.  A.,  Mathematical  Statistics  :  Basic  Ideas  and 
Selected  Topics,  Holden  -  Day  Inc.,  San  Francisco,  California,  1977. 

5.  Larson.  H.  J..  Introduction  to  Probability  Theory  and  Statistical  Inference,  3rd  ed., 
John  Wiley  and  Sons,  Inc.,  1982. 

6.  Duncan,  A.  J.,  Quality  Control  and  Industrial  Statistics,  5th  ed.,  Irwin,  Inc.,  1986. 

7.  Ouenouille,  M.  H.,  Notes  on  Bias  in  Estimation,  Biometrika,  43.  1956,  pp.  353  - 
360. 

8.  Tukey.  J.  W.,  Bias  and  Confidence  in  Not  Quite  Large  Samples,  Ann.  Math. 
Stat.,  29,  1958,  p.  614. 

9.  Miller,  R.  G.,  The  Jackknife  -  a  Review,  Biometrika,  6.  1974,  pp.  1  -  15. 

10.  Efron.  B..  Bootstrap  Methods  :  Another  Look  at  the  Jackknife,  Ann.  Statist.  7, 
1979,  pp.  1-26. 

11.  Mosteller.  F.  and  Tukev,  J.  W.,  Data  Analvsis  and  Regression  .  Addison- Weslev, 
Inc.,  1971  * 

12.  Efron,  B.    and  Tibshirani,  R.,  The  Bootstrap  ,  Statistical  Science,  1986,  pp.  54-76. 

13.  Naval  Postgraduate  School  Technical  Report,  NPS  55-81-005.  The  New  Naval 
Postgraduate  School  Random  Number  Package- LLRAN DOM  II.  bv  P.  A.  W. 
Lewis  and  L.  Uribe,  Monterev,  California,  1981. 


53 


INITIAL  DISTRIBUTION  LIST 


No.   Copies 


1.  Defense  Technical  Information  Center  2 
Cameron  Station 

Alexandria,  Virginia  22304-6145 

2.  Librarv.  Code  0142  2 
Naval  "Postgraduate  School 

Monterey,  California  93943-5002 

3.  Professor  P.AJacobs.  Code  55JC  2 
Department  of  Operations  Research 

Naval  Postgraduate  School 
Monterey.  California  93943-5000 

4.  Professor  Donald  Barr,  Code  55BN  1 
Department  of  Operations  Research 

Naval  Postgraduate  School 
Monterey,  California  93943-5000     • 

5.  Professor  Peter  Purdue.  Code  55  1 
Department  of  Operations  Researa. 

Naval  Postgraduate  School 
Monterey,  California  93943-5000 

6.  Librarv.  P.O.Box  77  2 
Pcstal'Code  130-09 

Gong,  Neung  Dong.  Dobong  Goo 
Seoul.  Republic  of  Korea 

7.  Major.  Kim.  Seung  Woong  15 
Postal  Code  151 

10-481  SinlimBon  Dong,  Gwanak.  Goo 
Seoul,  Republic  of  Korea 


54 


18i 


DUDLEY  KNOX  LIBRARY 
NAVAL  POSTGRADUATE  SCHOOL 
MONTEREY,  CALIFGi  943-8008 


Thesis 

1/A1R35  Kim 

i       A  simulation  of  esti- 
mates of  a  fi^t  passage 
time  distribution  for  a 
semi-Markov  process. 


